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INTRODUCTION 


Considerable  work  has  been  performed  over  the  past  15  years  on  the 
development  of  orthotropic  photoelastic  material  and  the  development  of 
stress-optic  laws  which  characterize  the  relations  between  the  optical 
and  mechanical  responses  of  the  material.  Some  plane  stress  elasticity 
problems  have  been  solved  with  the  material  but  these  problems  have  been 
more  or  less  of  a  demonstration  nature,  designed  to  show  the  extent  to 
which  orthotropic  photoelasticity  works.  In  many  cases  the  problems 
solved  were  such  that  the  load  was  applied  parallel  to  or  perpendicular 
to  the  fiber  direction,  the  principal  stress  system  coincided  with  the 
principal  material  system,  the  principal  stress  directions  were  known, 
the  stresses  were  computed  along  paths  parallel  or  perpendicular  to  the 
fiber  direction,  or  the  vast  majority  of  the  boundaries  were  traction 
free  so  that  2  out  of  the  3  stresses  were  known  with  certainty  over  much 
of  the  boundary.  Oftentimes  the  problems  were  chosen  because  there  ex¬ 
isted  known  solutions,  or  good  first  approximations,  to  the  state  of 
stress  in  the  body.  The  stresses  calculated  using  orthotropic  photo¬ 
elasticity  were  compared  to  the  known  solution.  These  types  of  exerci¬ 
ses  were  necessary  in  the  development  stages  of  this  new  technology. 
However,  to  have  the  same  utility  as  classical  isotropic  photoelastici- 
ty,  orthotropic  photoelasticity  must  be  applicable  to  a  wide  variety  of 
situations.  Ultimately,  orthotropic  photoelasticity  must  be  able  to  pre¬ 
dict  the  state  of  stress  when  the  applied  load  is  at  an  arbitrary  angle 
relative  to  the  fibers,  when  the  principal  stress  directions  are  vari¬ 
able  from  point  to  point  and  are  unknown,  and  when  the  stresses  along 
arbitrary  paths  are  required.  In  addition,  if  residual  birefringence 
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exists  in  the  material  due  to  the  manufacturing  process,  its  effect  on 
interpreting  photoelastic  data  must  be  taken  into  account. 

This  study  was  initiated  to  assess  the  ability  of  currently  accept¬ 
ed  orthotropic  stress-optic  laws  to  accurately  predict  the  state  of 
plane  stress  in  currently  available  orthotropic  photoelastic  material. 
The  ability  to  be  able  to  examine  a  variety  of  loading  directions,  rela¬ 
tive  to  the  axes  of  orthotropicity ,  was  of  prime  concern.  In  addition, 
it  was  important  to  be  able  to  compute  the  stresses  at  more  or  less  ar¬ 
bitrary  locations  in  the  interior  of  the  body.  In  classical  isotropic 
photoelasticity,  if  one  of  the  principal  stresses  is  known,  the  other 
principal  stress  can  be  computed  directly  from  the  isochromatic  fringe 
pattern.  This  is  always  the  case  when  examining  the  stresses  at  a  trac¬ 
tion-free  boundary.  For  this  case  the  shear  and  normal  stresses  are 
zero.  With  the  shear  stress  zero,  the  principal  stress  directions  are 
tangent  and  normal  to  the  traction-free  edge.  With  the  normal  stress 
zero,  the  only  unknown  is  the  principal  stress  tangent  to  the  free  edge. 
Using  the  stress-optic  coefficient  of  the  material  and  the  isochromatic 
fringe  count,  this  non-zero  principal  stress  can  be  computed.  In  the 
interior  of  a  body,  both  the  isochromatic  and  isoclinic  fringe  patterns 
are  used  to  compute  stresses.  Unless  a  third  point-by-point  experiment¬ 
al  measurement  is  used  in  conjunction  with  the  two  photoelastic  fringe 
measurements,  an  auxiliary  numerical  scheme  is  generally  used  to  compute 
the  complete  state  of  stress  in  the  interior  of  the  body.  This  auxili¬ 
ary  numerical  scheme  is  usually  in  the  form  of  the  numerical  integration 
of  the  plane-stress  equilibrium  equations  or  the  compatibility  equations 
along  a  path  from  a  boundary,  or  other  region  of  known  conditions,  to 
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the  interior  point  of  interest.  For  simplicity,  this  path  is  usually  a 
straight  line. 

This  report  begins  by  developing  the  stress-optic  law  for  orthotro¬ 
pic  materials.  The  laws  are  developed  on  a  physical  basis.  Some  of  the 
development  draws  from  the  past  work  in  orthotropic  photoelasticity  and 
can  be  considered  as  a  summary  and/or  a  reformating  of  results.  How¬ 
ever,  early  in  this  study  it  became  obvious  that  the  techniques  used  to 
fabricate  the  photoelastic  material  could  produce  a  residual  birefrin¬ 
gence.  With  no  applied  loads,  there  could  be  an  isochromatic  fringe 
count  of,  say,  0.1.  There  would,  of  course,  be  an  isoclinic  associated 
with  this  value.  The  development  which  follows  includes  the  effect  of 
residual  birefringence  in  the  stress-optic  law.  Conceptually,  account¬ 
ing  for  residual  birefringence  effects  is  simple.  Accounting  for  resid¬ 
ual  effects  operationally  is  also  quite  simple.  Since  birefringence  is 
associated  with  tensor  quantities,  keeping  track  of  residual  effects  is 
bookkeeping. 

To  have  the  ability  to  generate  a  wide  variety  of  stress  states  ex¬ 
perimentally  and  at  the  same  time  to  have  a  good  approximation  to  the 
state  of  stress,  an  unidirectional  disk  in  diametral  compression  was 
chosen  as  the  study  specimen.  With  such  a  specimen  the  fibers  could  be 
oriented  at  arbitrary  directions  relative  to  the  applied  diametral  load 
and  a  wide  variety  of  stress  states  could  be  produced.  In  addition,  the 
stress  state  along  arbitrary  diameters  or  chords  could  be  computed  using 
the  photoelastic  data  and  an  auxiliary  numerical  scheme.  It  was  decided 
to  refrain  from  using  a  third  point-by-point  experimental  measure  to 
compute  the  stresses.  The  most  accurate  point-by-point  measure  is  to 
use  interferometry  to  measure  the  thickness  change  of  the  material. 
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This  requires  optical  equipment  beyond  the  standard  polariscope  set-up. 
It  was  a  goal  to  keep  this  study  within  the  context  of  common  polari- 
scopssc  Oblique  incidence  techniques  were  not  considered  because  of  the 
poor  fringe  resolution  associated  with  orthotropic  photoelastic  materi- 

C  I  S  o 

The  disk  geometry  was  chosen  for  sever'al  other  reasons.  First j,  it 
■•■as  an  easy  specimen  geometry  to  rnake.  Second^  and  more  importantly^  it 
did  not  use  much  of  the  valuable  photoelastic  material  that  was  avail¬ 
able,  One  specimen  could  be  used  for  all  load  orientations.  Finally^ 
the  disk  geometry  is  an  ideal  specimen  for  calibration  of  orthotropic 
photoelastic  material.  In  theory,  according  to  the  state-of-the-art 
stress-optic  lav'S,  three  tests  must  be  made  on  an  orthotropic  photoelas¬ 
tic  material  to  completly  characterize  its  stress-optic  relations.  In 
the  past,  three  tensile  tests  have  been  conducted.  One  tensile  test  is 
conducted  on  a  0°  specimen,  one  test  on  a  90°  specimen,  and  one  test  on 
a  45°  specimen.  The  experimenter  can  judiciously  choose  the  length  and 
'Width  of  the  tensile  specimens  but  the  specimens  do  tend  to  consume  the 
valuable  photoelastic  material.  Since  the  material  has  some  variability 
in  its  mechanical  and  optical  properties,  it  is  a  good  idea  to  calibrate 
the  actual  panel  from  which  experimental  models  are  to  be  made.  Al¬ 
though  three  is  a  minimum,  conceivably  more  tensile  specimens,  say  at 
30°,  60°,  and  75°,  could  be  tested  and  the  calibration  could  be  done 
v/ith  some  form  of  least-squares  approach  to  the  data  reduction.  This 
would  be  done,  of  course,  at  the  expense  of  valuable  photoelastic  mate¬ 
rial,  This  has  not  been  done,  rather,  the  minimum  three  tests  have  gen¬ 
erally  been  the  standard.  With  a  single  disk,  theoretically,  an  infi¬ 
nite  number  of  stress  states  is  possible.  Practically  speaking,  with 
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various  disk  orientations  relative  to  the  load  direction,  and  with  pho¬ 
toelastic  data  taken  at  several  points  on  the  disk,  15-20  sets  of  stress 
and  corresponding  photoelastic  fringe  information  could  be  readily  ob¬ 
tained  and  a  least-squares  calibration  could  be  conducted  to  determine 
the  stress-optic  characterization  of  the  material. 

Experimentally  there  were  many  compelling  reasons  for  choosing  the 
disk  geometry.  In  addition,  the  theoretical  elasticity  solution  to  the 
disk  in  diametral  compression  could  be  derived  from  the  complex  variable 
approach.  The  complex  variable  solution  is  presented  in  an  appendix  and 
will  not  be  discussed  further  except  as  to  quote  stress  predictions  from 
the  solution. 

This  report  summarizes  the  work  performed  using  the  disk  as  a  study 
specimen.  The  work  followed  the  basic  line  of  reasoning:  If  a  particu¬ 
lar  stress-optic  law  for  the  orthotropic  photoelastic  material  is  true, 
and  if  the  calibration  of  the  material  is  done  using  three  tensile 
tests,  then  two  types  of  experiments  can  be  done  to  test  the  validity  of 
the  stress-optic  law.  Using  the  calibration  information,  the  stress- 
optic  law,  and  the  theoretical  elasticity  solution,  the  isochromatic  and 
isoclinic  fringe  patterns  in  the  disk  can  be  predicted.  By  examining 
different  load  orientations,  relative  to  the  fiber  directions,  a  wide 
variety  of  spatial  variations  of  fringe  patterns  can  be  generated.  This 
wide  variety  of  predicted  fringe  patterns  can  then  be  compared  to  the 
experimentally  observed  fringe  patterns.  This  is  a  legitimate  test  of 
the  assumed  stress-optic  behavior  and  the  calibration.  These  so-called 
"forward  tests"  were  conducted  and  are  reported.  Real  interest,  of 
course,  is  in  starting  with  a  specific  fringe  pattern  and  working  back¬ 
wards  to  obtain  the  stress  state.  This  is  the  essence  of  photoelastici- 
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ty »  A  numerical  scheme  was  previously  developed  for  studying  the  stres¬ 
ses  in  the  inner  laps  of  double-lappeds  double-pin  connector  models  con¬ 
structed  from  acrylic  and  isotropic  photoelastic  material.  This  scheme 
was  based  on  the  stress-optic  law  for  classic  isotropic  materials  and  a 
finite-difference  form  of  the  plane-stress  equilibrium  equations.  The 
equilibrium  equations  were  chosen  as  the  auxiliary  conditions  to  obtain 
the  three  components  of  stress  at  interior  points  in  the  model.  This 
scheme  was  modified  to  accomodate  an  orthotropic  stress-optic  law.  The 
stress-optic  calibration  information  from  auxiliary  calibration  tests, 
tne  isochromatic  and  isoclinic  fringe  count  at  the  various  mesh  points 
of  the  finite-difference  grid,  and  the  finite-difference  equilibrium 
equations  were  used  to  experimentally  determine  the  stresses  in  the 
disk.  These  stresses  were  compared  with  the  stresses  computed  from  the 
complex  variable  elasticity  solution.  These  stress  prediction  tests 
could  be  considered  as  "backwards  tests,"  To  put  the  results  into  con¬ 
text,  numerical  and  experimental  results  for  disks  made  from  isotropic 
photoelastic  material  are  compared.  Also  numerical  studies  were  conduc¬ 
ted  to  determine  the  effect  of  isocline  measurement  error  on  the  experi¬ 
mental  results. 

Throughout  the  report  specific  references  to  previous  investiga¬ 
tions  are  made.  These  are  noted  in  a  list  of  references.  However,  for 
completeness,  a  bibliography  is  included  at  the  end  of  this  report  to 
mention  those  investigations  not  specifically  cited. 
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ORTHOTROPIC  PHOTOELASTICITY 


Bi refringence 

In  many  transparent  solids,  the  Index  of  refraction*  can  be  consid¬ 
ered  as  a  quantity  which  is  independent  of  the  direction  of  propagation 
and  the  plane  of  polarization  of  the  transmitted  light.  However,  in 
some  transparent  solids  the  speed  of  propagation  of  light  depends  both 
on  the  direction  of  propagation  and  the  plane  of  polarization  of  the 
light.  Within  such  a  solid  there  will  be  a  direction  of  propagation  and 
a  plane  of  polarization  for  which  the  light  travels  the  fastest.  There 
is  also  a  direction  of  propagation  and  plane  of  polarization  for  which 
light  travels  the  slowest.  Because  of  the  orthogonality  and  symmetry  of 
the  molecular  structure  and  crystals  making  up  these  solids,  many  of 
these  directions  and  planes  are  perpendicular  to  each  other.  This  idea 
of  multiple  propagation  speeds  and  directional -dependence  of  the  speed 
of  propagation  is  termed  bi refri ngence.  In  the  discussions  here,  the 
birefringence  phenomenon  will  be  considered  something  that  can  occur  in 
an  unstressed  material  and  also  something  that  can  be  induced  in  a  ma¬ 
terial  by  application  of  loads.  Because  of  the  orthogonality  of  the 
directions  associated  with  the  different  propagation  speeds,  the  indices 
of  refraction  can  be  associated  with  quantities  which  can  be  regarded  as 
tenso'^s.  Tensor  quantities  have  principal  values  and  the  directions  as¬ 
sociated  with  these  values  are  quite  often  perpendicular  to  each  other. 
The  ideas  of  extremums  in  the  propagation  speeds,  tensors,  planes  of  po¬ 
larization,  etc.,  can  all  be  used  to  put  the  ideas  of  classical 


*index  of  refraction  of  a  medium,  n,  is, defined  as  the  speed  of  light  in 
a  vacuum,  Cq,  divided  by  the  speed  of  light  in  the  medium,  c. 
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isotropic  photoelasticity  on  a  formal  and  rigorous  foundationo  This  has 
been  done  in  various  texts  and  papers.  However^  if  one  is  trying  to  de- 
terminej  for  examples  the  stress  concentration  factor  at  the  base  of  a 
fillet  in  a  model  representing  a  metal  components  the  formalism  is  not 
required.  All  that  is  required  is  to  know  that  the  isochromatic  fringe 
number  times  some  a  priori  determined  numerical  factor  gives  the  value 
of  the  stress  at  the  base  of  the  fillet.  However s  if  one  is  dealing 
with  a  new  photoelastic  phenomenon  such  as  orthotropic  photoel asticitys 
the  formalism  helps  guide  the  interpretation  of  the  photoelastic  data. 
The  formalism  also  provides  insight  into  the  potential  of  the  new  phe¬ 
nomenon.  In  this  reports  some  of  the  formalism  will  be  used  to  estab¬ 
lish  ideas  regarding  orthotropic  photoelasticity.  The  dielectric  tensor 
and  its  relation  to  the  indices  of  refraction  will  be  introduced.  Then 
the  notion  of  the  birefringence  tensor  will  be  discussed.  The  introduc¬ 
tion  of  the  birefringence  tensor  is  a  convenient  abstraction  and  is  not 
at  all  necessary  in  order  to  understand  the  implications  of  orthotropic 
photoelasticity.  However,  many  past  investigators  have  used  the  idea  of 
a  birefringence  tensor  and  so  it  will  be  used  here  to  illustrate  how  it 
fits  into  the  nomenclature  and  mechanics  of  orthotropic  photoelasticity. 
These  ideas  will  be  used  to  derive  relations  betvveen  the  applied  stress 
and  the  observed  optical  response.  These  relations  will  be  such  that 
they  can  be  used  in  a  numerical  scheme  to  compute  the  stresses  in  a 
loaded  photoelastic  model.  To  begin,  however,  the  approach  used  by  many 
texts  to  start  the  discussion  of  isotropic  photoelasticity  will  be  re¬ 
produced.  The  expression  for  the  relative  retardation  will  be  derived. 
This  approach  will  introduce  the  relations  between  the  wavelength  of 
light,  the  model  thickness,  the  indices  of  refraction,  and  the  slowing 
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of  the  light  passing  through  photoelastic  material.  For  those  not 
familiar  with  even  isotropic  photoelasticity,  this  approach  will  be 
quite  palatable  because  it  is  based  on  the  ideas  of  travel  times  of  a 
light  wave,  optical  path  lengths,  etc.  Also,  since  many  texts  use  this 
idea,  it  will  be  a  point  of  reference  if  one  is  consulting  these  texts. 
For  those  familiar  with  isotropic  photoelasticity,  the  derivation  will 
be  a  review  and  it  will  provide  some  commonality  between  isotropic  and 
orthotropic  photoelasticity.  The  relative  retardation  is  introduced  be 
cause  it  is  one  of  the  two  basic  quantities  measured  directly  by  a 
standard  polariscope.  The  other  quantity,  which  will  be  introduced 
after  the  relative  retardation  is  discussed,  is  the  isocline. 

Relative  Retardation 

Consider  a  transparent  solid  in  the  form  of  a  thin  plate,  of  thick 
ness  h,  oriented  perpendicular  to  the  direction  of  propagation  of  a 
light  vector.  As  shown  in  fig.  1,  a  light  vector  is  impinging  on  or  in 
cident  to  the  plate.  If  the  light  vector  is  polarized,  it  can  be 
thought  of  as  oscillating  harmonically  in  a  single  plane  as  it  propa¬ 
gates  toward  the  plate,  as  shown  in  fig.  2,  If  the  plate  is  made  of  a 
material  with  indices  of  refraction  which  are  dependent  on  the  plane  of 
polarization  of  the  incident  light,  then  it  is  possible  to  identify  a 
plane  of  polarization  in  which  the  light  travels  through  the  plate  with 
the  fastest  propagation  velocity.  It  is  also  possible  to  identify  a 
plane  of  polarization  in  which  the  light  passes  through  the  plate  with 
the  slowest  propagation  velocity.  These  directions,  as  shown  in  fig,  3 
are  referred  to  as  fast  and  slow  directions,  respectively.  The  descrip 
tions  fast  and  slow  can  be  misleading.  In  a  transparent  solid  light 
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actually  slows,  when  compared  to  its  speed  in  air,  independent  of  the 
plane  of  polarization.  However,  in  the  so-called  fast  direction  the 
light  slows  the  least. 

The  polarized  light  that  is  impingent  on  the  plate  in  some  direc¬ 
tion  other  than  the  fast  or  slow  direction,  as  shown  in  figure  4,  can  be 
decomposed  into  vector  components  in  the  fast  direction  and  slow  direc¬ 
tion,  respectively.  Referring  to  in  fig.  5,  due  to  differences  in  the 
speed  of  propagation  in  the  two  directions,  these  two  components  are 
out-of-phase  with  each  other  when  they  emerge  from  the  plate.  Since  the 
components  are  out  of  phase  they  can  be  made  to  interfere  constructively 
and  destructively.  The  specific  conditions  required  to  obtain  inter¬ 
ference  is  a  function  of  the  optical  apparatus  used  for  viewing  the  phe¬ 
nomenon.  The  various  optical  devices,  or  polariscope  arrangements,  and 
the  specific  relations  for  interference  of  these  out-of-phase  components 
are  discussed  in  [1].  It  is  this  phase  difference  which  is  the  basis 
for  using  the  birefringent  phenomenon  to  an  advantage.  An  expression 
for  this  phase  difference  is  developed  in  the  next  section. 

Phase  Difference 

If  the  index  of  refraction  in  the  fast  direction  is  denoted  as  n2, 

then 

where  Cq  is  the  speed  of  light  in  a  vacuum  and  C2  is  the  speed  of  light 
in  the  plate  if  the  polarized  light  vector  is  aligned  with  the  fast,  2, 
direction.  The  index  of  refraction  in  the  slow,  1,  direction  is  simi¬ 
larly  defined,  i.e. 
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n 


1 


The  hats  over  the  indices  of  refraction  signifies  that  these  are  the 
extremum  values  for  all  orientations  of  the  plane  of  polarization.  This 
hat  notation  will  soon  be  used  to  denote  the  principal  values  of  a  ten¬ 
sor  associated  with  the  transmission  of  light  through  the  plate.  The 
principal  values  of  this  tensor  and  the  extremums  in  the  values  of  the 
indices  of  refraction  are  related. 

For  light  propagating  with  the  fast  speed,  to  pass  through  thick¬ 
ness  h,  the  elapsed  time,  t2,  is 


Similarly,  for  light  propagating  with  the  slow  speed,  the  elapsed  time, 
tj,  to  pass  through  thickness  h  is 


If  a  vacuum  occupied  that  region  of  space  instead  of  a  birefringent 
medium,  then  the  elapsed  time  to  move  through  the  distance  h  would  be 


t 


0 


(5) 


Figure  6  illustrates  this  thickness  geometry.  The  time  difference  be¬ 
tween  when  the  fast  light  emerges  from  the  plate  and  when  the  slow  light 
emerges  from  the  plate.  At,  is  given  by 


16 


At 


(6) 


Using  eqSo  1-4  in  eq„  6  leads  to 

At  “  —  (Oh  —  Op)  o 
0 


(7) 


In  terms  of  distance  (e.g.  meters s  inches j  etc„)s  the  two  light  vectors 
are  out  of  phase  by  an  amount 

c^At  =  h  (n^  -  n^)  .  (8) 


Relative  to  the  wavelength  of  light  g  x,  entering  the  plate^  the  phase 
difference  5  in  cycles  $  between  the  two  light  vectors  is 


CoAt 


(9) 


To  put 
by  ?iig 


this  phase  difference  in  terms  of  radians,  multiply  eq, 
i  .00 

_  2TrcQAt  _  , 

^  X  “  X  ~  ^^2^  ° 


(9) 

(10) 


Thus  the  phase  difference,  a,  is  expressed  in  terms  of  the  wavelength  of 
incident  light,  the  plate  thickness,  and  the  difference  in  the  principal 
indices  of  refraction.  The  quantity  A  is  the  relative  retardation  be- 
tween  the  two  components  of  reemerging  polarized  light. 

In  this  discussion  the  quantity  n^  is  arbitrarily  associated  with 
the  direction  that  transmitts  the  light  the  fastest  and  n^^  is  associated 
with  the  direction  that  transmitts  the  light  the  slowest.  Thus 
(ni  -  n^)  is  a  positive  quantity. 
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Causes  of  Bi'refringence  in  Solids 

The  interaction  of  the  electromagnetic  field  of  the  light  with  the 
electromagnetic  fields  associated  with  the  bonding  of  the  atoms  in  a 
transparent  solid  changes  the  speed  of  light  from  the  vacuum  condition. 
In  addition,  anything  that  changes  the  spatial  distribution  and/or 
strength  of  the  bonding  electric  fields  in  a  solid  will  affect  the  pro¬ 
pagation  speed  of  light  within  the  solid.  A  way  to  change  the  bonding 
electric  fields  in  a  solid  is  to  physically  move  the  molecules  of  the 
solid  relative  to  each  other.  Such  a  condition  could  be  caused  by  ap¬ 
plying  a  mechanical  load  (stress)  to  the  solid.  If  the  solid  is  elas¬ 
tically  isotropic,  and  if  there  is  a  stress  applied  in  a  particular 
direction,  then  the  strength  and  distribution  of  the  bonding  fields 
would  be  different  in  different  directions.  In  this  case  the  electric 
field  of  the  incident  light  would  be  acted  upon  by  one  field  strength  if 
it  was  polarized  in  one  direction  and  it  would  be  acted  upon  by  another 
field  strength  if  it  was  polarized  in  another  direction.  This  leads  to 
di recti onally-dependent  propagation  speeds.  This  phenomenon  would  also 
occur  if  the  material  was  elastically  orthotropic.  In  fact,  even  with 
the  solid  in  its  natural  (unstrained  or  unstressed)  state,  an  orthotro¬ 
pic  solid  has  different  bonding  field  strengths  in  different  directions. 
Thus  it  is  possible  to  conceive  of  multiple  indices  of  refraction  in  an 
orthotropic  media  with  the  molecules  in  their  natural  (unstressed)  posi¬ 
tions.  In  the  context  of  these  discussions,  the  multiple  indices  of  re¬ 
fraction  with  no  load  applied  can  be  considered  a  residual  difference  in 
the  indices  of  refraction  or  residual  birefringence.  If  a  mechanically 
induced  alteration  of  the  bonds  occurs,  and  if  superposition  of  bire- 
fringent  effects  is  assumed  valid,  then  the  total  differences  in  indices 
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rion  would  be  the  sum  of  the  effects  of  the  residual  bire- 
and  the  effects  of  the  niechanically  induced  bi refringence:, 


1 1s  of  Birefringence  in  Solids 

Tiie  fundamental  quantities  for  the  study  of  bi  ref  ri  ngsnce  in  solids 


Ui* 


ro 


dielectric  constants  of  the  material.  The  plural  connotation  in 
d  'constants'  is  important  because  oftentimes  the  dielectric 
las  of  a  material  are  thought  of  as  being  described  by  a  single 


:  In  general j  the  dielectric  properties  of  a  material  are  dire 
ly-dvspandent  and  they  can  be  represented  as  a  tensor  called  the 
The  dielectric  tensor  relates  the  electric  field 
ti;  vector  to  the  electric  displacement  vector.  The  index  of  re 
0,1  of  a  solid  is  related  to  the  dielectric  tensor.  Due  to  the 


di  - 


rcl ly-dependent  nature  of  a  tensor,  the  speed  of  propagation  of 
v:i  a  transparent  solid  depends  strongly  on  the  di  recti  onal -depen- 
of  the  dielectric  properties  of  the  solid.  For  the  case  of  a  thin 
.'•aferred  to  here  by  analogy  to  'plane  stress'  as  'plane  birefrln- 
■  tlie  dielectric  tensor  can  be  written  as 


thie  above  notation,  the  dielectric  properties  are  being  referred  to 
y  coordinate  system.  The  x-y  coordinate  system  is  a  global  coordi- 
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nate  system  with  the  x  and  y  axes  in  the  plane  of  the  plate.  At  the  mo¬ 
ment  the  orientation  of  the  x  and  y  axes  within  the  plane  are  not  impor¬ 
tant,  Later,  stresses  will  be  computed  in  the  x-y  system  and  its  orien¬ 
tation  will  be  specified.  Figure  7  depicts  the  coordinate  system.  The 
light  can  be  considered  to  be  propagating  into  the  page  and  perpendicu¬ 
lar  to  it. 

In  an  electrically  isotropic  material  the  dielectric  tensor  takes 
the  form 

e  ) 

e[  .  (12) 

0  \ 

A  vacuum,  or  free  space,  is  isotropic  and  so  the  dielectric  tensor  is 
given  by 


the  subscript  'o'  denoting  free  vacuum.  The  fast  and  slow  indices  of 
refraction  of  a  solid  are  related  to  the  principal  values  of  the 
dielectric  tensor.  The  principal  values  of  the  dielectric  tensor  are 
given  by 

(14) 


The  fast  and  slow  indices  of  refraction  are  defined  [2]  to  be 


For  the  transparent  solids  considered  here  and  for  the  levels  of 
birefringence  encountered  in  photoelasticity,  certain  simplifying  ap¬ 
proximations  can  be  made.  If  e^^g  is  defined  as 


then  eq.  14  becomes 


where 


(17) 

(18) 


It  is  shown  in  Appendix  A  that  Ae  is  an  order  of  magnitude  smaller  than 
®ave*  Then,  using  the  binomial  approximation,  eqs.  15  and  16  become 


The  relative  retardation  can  be  related  to  the  dielectric  constants 
by  substituting  eqs.  20  and  21  into  eq.  10.  This  relation  is  important 
because  the  application  of  stress  directly  affects  the  dielectric  con¬ 
stants  of  a  transparent  material.  The  above  substitution  results  in 
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A 


(22) 
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where  is  the  average  index  of  refraction  of  the  medium  and  is 
defined  as 


It  should  be  mentioned  at  this  point  that  it  has  been  assumed  that 
the  magnetic  permeability,  y,  has  been  assumed  to  be  an  isotropic  prop¬ 
erty  of  the  material.  Furthermore,  the  permeability  is  assumed  to  be 
equal  to  the  permeability  in  a  vacuum,  Generally  the  velocity  of 
propagation  of  electromagnetic  waves  depends  on  the  permeability  as  well 
as  the  dielectric  constant.  For  an  isotropic  material,  the  velocity  of 
propagation  is  proportional  to  l/ZeiT  .  In  a  vacuum  the  velocity  becomes 
index  of  refraction  of  a  medium  is  thus  given  by  thus 


If  y  is  assumed  to  be  equal  to  y^,  then  the  index  of  refraction  is 
related  only  by  the  ratios  of  the  dielectric  constants. 

For  an  orthotropic  material,  particularly  fiber  reinforced  compos¬ 
ite  materials,  there  is  a  coordinate  associated  with'  the  stiff  (fiber) 
direction  and  a  coordinate,  perpendicular  to  the  first,  associated  with 
the  soft  (matrix)  direction.  Much  of  the  analysis  of  fiber-reinforced 
composites  uses  these  directions  as  the  reference  for  defining  stress, 
strains,  elastic  constants,  failure  stresses,  etc.  In  orthotropic  pho¬ 
toelasticity  it  is  also  convenient  to  use  this  coordinate  system  as  a 
basis.  Herein,  tnese  directions  will  be  referred  to  as  IM  and  2M„ 
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These  directions  will  refer  to  the  fiber  and  matrix  direction,  respec¬ 
tively,  with  the  M  denoting  that  it  is  a  material  coordinate  system. 

When  computing  stresses  in  the  x-y  system,  the  angle  the  +1M  axis  makes 
with  to  the  +x  axis  will  be  denoted  This  coordinate  system  orienta¬ 
tion  is  shown  in  fig.  8.  In  the  material  coordinate  system,  the  dielec¬ 
tric  tensor  can  be  identified  as 


Birefringence  in  Fiber-Reinforced  Solids 

Because  of  differences  in  the  thermoelastic  properties  of  the  fiber 
and  matrix  constituents  in  fiber-reinforced  composites,  the  fabrication 
of  fiber-reinforced  composites  can  lead  to  residual  stresses  in  the  fi¬ 
bers  and  the  matrix.  Since  the  dielectric  properties  of  solids  can  de¬ 
pend  on  the  stress  state  in  the  solid,  these  residual  stresses  can  lead 
to  residual  birefringence.  As  mentioned  in  the  Introduction,  such  bire¬ 
fringence  is  considered  here.  In  the  material  coordinate  system,  resid¬ 
ual  dielectric  effects  generally  have  three  components. 


The  superscript  R  identifies  properties  associated  with  residual  ef¬ 
fects.  There  is  no  a  priori  reason  to  believe,  in  general,  that  the 
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principal  directions  of  the  residual  dielectric  tensor  are  aligned  with 
the  material  directions.  There  is,  however,  a  principal  direction  asso 
ciated  with  the  residual  tensor,  i.e.  directions  associated  with  the 
fastest  and  slowest  propagation  directions  for  the  solid  when  there  is 
no  mechanical  load  applied.  Assume  that  the  principal  residual  direc¬ 
tions  are  rotated  from  the  material  directions  by  an  amount  9^.  The  or 
ientation  of  the  1R-2R  system  relative  to  the  1M-2M  system  is  shown  in 
fig.  9.  The  residual  dielectric  tensor  principal  coordinate  system  is 
identified  by  IR  and  2R.  The  principal  values  of  the  residual  dielec¬ 
tric  tensor  are  given  by  e^j^  and  e2|^.  The  representation  of  the  residu 
al  dielectric  tensor  in  the  principal  residual  coordinate  system  is 


The  components  of  the  residual  dielectric  tensor  in  the  material 
system,  eq.  26,  can  be  related  to  the  principal  values  of  the  residual 
dielectric  tensor,  eq.  27,  by  the  usual  two-dimensional  tensor  trans¬ 
formations,  i.e. 
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The  minus  sign  is  left  with  the  angle  as  a  reminder  that  in  these  equa¬ 
tions  one  is  transforming  back  to  the  material  system.  If  the  compo¬ 
nents  of  the  residual  dielectric  tensor  are  known  in  the  material  sys¬ 
tem,  the  principal  values  and  principal  directions  of  the  residual  diel¬ 
ectric  tensor  can  be  determined  by  the  familiar  formulas 


26^ 

tan  (2ej^)  =  •  (2®)® 

®1M  "  ®2M 

If  loads  are  applied  to  the  solid  and  there  are  residual  birefrin- 
gent  effects  present,  then  there  will  exist  a  state  of  total  birefrin- 
gence.  This  situation  is  associated  with  the  total  dielectric  tensor. 
The  total  dielectric  tensor  will  have  principal  values  and  principal 
directions.  The  coordinate  system  of  the  principal  total  dielectric 
tensor  will  be  denoted  as  IT  and  2T,  the  T  signifying  association  with 
total  effects.  The  principal  values  of  the  total  dielectric  tensor  can 
be  written  as  , 

I  ^ 

yn 

|®2T 
( ° 

In  the  material  coordinate  system,  the  components  of  the  total  dielec¬ 
tric  tensor  are  given  by 
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(30) 


eT 
^2M 

eT 
12M 

Ths  angle  betv/een  the  material  coordinates  and  the  coordinate  system  of 
the  principal  total  dielectric  tensor  is  denoted  as  e-j.  and  is  illustra¬ 
ted  in  f1go  lOo  The  components  of  the  total  dielectric  tensor  in  the 
material  coordinate  system^  eq„  30,  are  related  to  the  principal  values 
of  the  total  dielectric  tensor,  eq,  29,  by 


T 

.IT  ■  ®2T^ 

.  2  J 

cos  (-26^) 

(31)a 

^2M 

r®lT  ®2T^ 

^  2  J  ■ 

.®1T  ■  ®2Tv 

1  2  J 

cos  (-26-j-)„ 

(31)b 

e> 

12M 

sin  (-2ej) 

(31)c 

Conversely,  knowing  the  components  of  the  total  dielectric  tensor  in  the 
material  system,  the  principal  values  and  principal  directions  of  total 
effects  can  be  computed  by 


®1T’'^2T 


e^  + 

•  IM  2M 


■)±  V( 


e^  + 


(e-i 


12M^ 


(31)d 


2e 


tan  (2©^)  =  —j- 


12M 


(e 


J 


(31)e 


IM 


rinally,  there  are  principal  values,  principal  directions,  and  ma¬ 
terial  system  components  of  the  dielectric  tensor  associated  with  the 
application  of  mechanical  loadSo  In  the  material  system,  the  components 

of  the  dielectric  tensor  due  to  the  applied  loads  are  given  by 

30 


Fig.  10  Orientation  of  1T-2T  syst 
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(32) 


The  mechanicany  induced  dielectric  tensor  has  its  own  principal  coordi¬ 
nate  system  la  -  2a„  The  principal  values  of  the  dielectric  tensor  due 
to  applied  loads  can  be  referred  to  this  system  and  can  be  written  as 


Assume  the  lo  -  2a  coordinate  system  is  oriented  at  an  angle  relative 
to  the  material  system„  Figure  11  shows  the  angle  9^  and  four  of  the 
five  coordinate  systems  discussedo  (The  x-y  coordinate  system  is  not  de¬ 
picted.)  Equations  similar  to  eqs.  28a  -  28e  and  31a~31e  can  be  written 
to  relate  ej^,  e^^,  and  e^.  These  equations  are 


e.  +  e^ 

■^IM  =  (-IT--)  "  "-OS  (-29J 
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;34)a 


e.  +  e^  e,  ~  e^ 


2 


(34)b 


"la  -  "2a 


'12M 


]  sin  (-29  ) 
J  ^  a 


'34)c 


and 
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Hypothesis  Regarding  Dielectric  Effects  Due  to  Applied  Lqadi 

It  will  be  hypothesised  here  that  the  dielectric  effects  due  to 
applied  loads  can  be  superimposed  onto  the  residual  dielectric  ef- 
fectSo  Due  to  the  tensorial  nature  of  dielectric  effects,  the  total 
dielectric  tensor  can  be  represented,  in  the  material  system,  as 


follows: 

\  / 

e'^  1 
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It  has  been  hypothesized  by  previous  investigators  [3,4]  that  the 
dielectric  tensor  due  to  the  applied  loads  is  related  to  the  applied 
stress  by  the  following  orthotropic  relation: 


(36) 


The  q-jj  are  constants  for  a  given  material  and  can  be  thought  of  as  the 
material's  stress-dielectric  coefficientSo  Equation  35,  representing 
the  total  dielectric  tensor  in  the  material  system,  can  be  written  as 
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Several  approaches  can  be  taken  to  check  the  validity  of  eq.  37  or 
to  check  the  consequences  of  it,  assuming  it  is  valid.  It  is  on  the 
basis  of  eq.  36  that  much  of  the  work  in  orthotropic  photoelasticity  has 
been  done  to  date.  The  hypothesis  represented  by  eq.  36  has  been  termed 
ad  hoc  [5]  because  there  is  no  a  priori  reason  to  believe  that  orthotro¬ 
pic  elastic  behavior  implies  orthotropic  optical  behavior.  The  crux  of 
the  discussion  about  the  validity  of  eq.  36  is  whether  or  not  it  corre¬ 
lates  with  experimental  observation.  Equation  35  has  not  been  verified 
to  any  great  extent.  Previous  investigators  [6, 7, 8, 9]  have  mentioned 
residual  birefringent  effects.  However  eq.  35,  and  the  just-derived 
relations  involving  the  residual  dielectric  tensor,  treat  residual 
effects  with  more  formalism. 


Consequences  of  the  Hypothesized  Relation 

Since,  as  stated  previously,  standard  polariscopes  measure  relative 
retardation.  A,  it  is  useful  to  relate  the  retardation  a  to  the  applied 
stresses  via  the  just  stated  hypotheses.  This  is  done  by  using  eq.  37 
in  eq.  22,  the  expression  for  a.  In  the  present  context  the  dielectric 
properties  in  eq.  22  refer  to  the  total  dielectric  properties,  i.e. 
residual  effects  plus  mechanically  induced  effects.  Thus  eq.  22  becomes 
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Xn 


2TTh 
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ave 


(38) 
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Since  1s  related  to  the  trace  of  the  total  dielectric  tensor ^  e^^^ 

is  independent  of  coordinate  system.  To  use  eq.  37  in  eq.  39 ^  it  is 
useful  to  write  the  residual  dielectric  effects  in  eq.  37  in  terms  of 
principal  residual  values.  This  can  be  done  using  eqs.  28d  -  23c, 

Doing  this^  eq,  37  becomes 
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Substituting  eq,  40  into  eq,  39  and  then  substituting  those  results  into 
eoo  38  leads  to  an  expression  for  A,  namely 


Xn  G  ^ 
ave  0 


\/{{qil  -  “  ^^22  ““  ^12^^2  “  ^2R^  cos(2ej^)}^ 


+  {2qccT^-jo  (Sin  -  SpR^  sin(29p)}2 


With  a  darkfield  polariscope  arrangement,  for  example,  fringes  occur 
when 

■|-  =  0 ,  IT ,  2ir ,  3ti  , . , , ,  Nir  ( 42 ) 

N  an  integer. 

Combining  eqs.  41  and  42  gives  the  condition  for  integer  fringe  orders 
in  a  darkfield  polariscope  due  to  applied  loads  and  residual  effects. 
Specifically,  integer  order  fringes  occur  when 


Nj  being  the  integer  fringe  order.  The  subscript  T  has  been  added  to  N 
to  signify  it  is  the  fringe  order  associated  with  total  birefringence, 
residual  and  mechanically  induced.  Practically  speaking,  Nj  does  not 
have  to  be  an  integer.  Standard  polariscopes  are  capable  of  displaying 
any  fringe  order,  fractional  or  integer.  In  any  case,  the  residual  ef¬ 
fects  and  the  applied  load  effects  combine  in  the  manner  as  indicated  by 
the  left  hand  side  of  eq.  43  to  produce  the  observed  fringe  order.  This 
is  a  direct  consequence  of  eq.  37  and  the  optics  of  a  darkfield  polari¬ 
scope.  Knowing  the  material  constants  q^-j  and  the  residual  effects,  eq. 
43  can  be  used  to  determine  the  fringe  order  produced  by  a  certain 
stress  state.  It  is  important  to  note  that  in  eq.  43  the  five  optical 
constants  q^-j  appear  as  differences  in  two  of  the  coefficients  of 
stress.  In  reality,  only  (q^  -  q2i),  (qi2  -  ^zz)*  and  q55  are  neces- 
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sary  for  determining  the  effects  of  applied  stress  on  fringe  order. 

This  constitutes  three  quantities  instead  of  five.  In  addition,  only 
the  difference  in  principal  values  of  the  residual  dielectric  tensor  is 
involved,  along  with  the  principal  residual  directions. 

It  is  legitimate  at  this  point  to  be  concerned  with  the  fact  that 
to  use  eq,  43,  such  unfamiliar  physical  constants  as  the  average  total 
index  of  refraction,  n"*"  ,  the  stress-dielectric  coefficients,  q^^,  the 

residual  dielectric  effects,  etc,  have  to  be  determined.  In  the  next 
section  many  of  these  constants  will  be  lumped  together  into  stress- 
optic  coefficients  and  another  tensor  called  the  birefringent  tensor. 
Then  instead  of  having  to  be  concerned  with  all  of  the  quantities  in  eq, 
43,  concern  will  be  with,  for  example,  the  slopes  of  the  fringe  order 
vs,  applied  stress  relations. 


Birefringence  Tensor 


To  examine  other  consequences  of  eq,  37,  it  is  convenient  to  define 


several  quantities  which  can  be  thought  of  as  components  of  another 
tensor.  These  components  combine  tensorially  to  produce  the  fringe  or¬ 


der  Ny.  Define 
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where  a^j 


h 


i,  j  =  1,  2,  6. 


(46) 


The  various  quantities  N  have  been  called  components  of  the  birefringent 
tensor  (as  opposed  to  components  of  the  dielectric  tensor).  Having  par¬ 
ticular  values  of  the  components  and  is  having  a  particu¬ 
lar  ^t^te_of_bij2efjri^^  Using  eqs.  46  and  47,  eq.  43  becomes 


This  equation  indicates  that  the  observed  fringe  count  Nj  is  due  in  part 
to  the  effect  of  normal  stresses,  and  N^j^,  the  shear  stress,  N®2|v|» 
and  the  residual  birefringence.  The  equation  indicates  that  the  fringe 
productions  due  to  normal  stresses,  and  N^j^,  subtract  algebraically 
and  the  difference  is  added  vectorially  to  the  fringe  production  due  to 
the  shear  stress,  2NJ2hi|.  Residual  effect  fringe  production  adds  alge¬ 
braically  to  the  fringe  production  due  to  applied  load.  Equation  49 
states  that  all  effects  contribute  to  the  observed  fringe.  In  standard 
polariscopes,  the  individual  components,  i.e.  N^j^,  N2|„],...etc.,  cannot 
be  observed.  Only  the  total  effect,  Nj,  is  observable.  In  a  sense. 


then,  the  quantities  N2[(^,...etc„  are  simply  convenient  abstractions 

which  essentially  replace  the  dielectric  tensor. 

Mohr's  Circle  of  Birefringence 

Equation  49  has  an  interesting  physical  interpretation.  Consider, 
for  the  moment,  the  case  where  the  residual  effects  are  zero.  Equation 
49  becomes 


If  a  Mohr's  circle  is  constructed  with  and  on  the  horizontal  ax¬ 
is  (i.e.5  where  a  and  a  are  located  on  Mohr's  circle  of  stress)  and 
X  y 

on  the  vertical  axis,  then  the  state  of  birefringence  due  to  the 
applied  load  can  be  located  on  these  axes.  The  state  of  birefringence 
is  shown  as  closed  circles  on  fig.  12a.  By  geometry,  the  quantity 

/  -  N® 

\/(— +  (Num)" 

is  the  radius  of  the  circle.  Then  twice  the  quantity  in  eq.  51  is  the 
diameter  of  the  circle.  Thus,  by  eq.  50,  with  no  residual  effects  pres¬ 
ent,  the  diameter  of  Mohr's  circle  of  birefringence  is  the  fringe  order. 
My,  seen  in  a  polariscope.  Equation  50  is  the  form  originally  intro¬ 
duced  by  Sampson  [10]  for  utilizing  the  idea  of  a  birefringent  tensor. 

He  essentially  postulated  eq.  50  and  worked  from  there.  Figure  12b 
shows  other  important  aspects  of  Mohr's  circle  of  birefringence. 

A  Mohr's  circle  interpretation  can  be  implied  even  in  the  presence 
of  residual  birefringence.  This  circle  is  Mohr's  circle  of  total  hire- 
fringence.  Referring  to  eq,  40,  define 


40 


T  T 

Construct  Mohr's  circle  by  locating  and  on  the  horizontal  axis 
T 

and  the  vertical  axis.  The  diameter  of  the  circle  is 


(53) 


Using  the  definitions  of  eqs.  45-47  and  substituting  into  eq.  53,  the 
diameter  expression  becomes 

-  Nji^)  +  (Njp  -  N2r)  cos(2eu))2  + 

{2N»2|,  +  (Ni„  -  Njr)  s1n(2e„)}2  . 

Thus  an  interpretation  of  eq.  49  is  that  the  diameter  of  Mohr's  circle 
of  total  birefringence  is  equal  to  the  total  fringe  order,  Nj,  observed 
in  the  polariscope.  Mohr's  circle  of  total  bi ref ri ngence  is  illustrated 
In  fig.  13,  Because  of  the  connection  with  Mohr's  circle,  it  is  obvious 
that  the  observed  total  fringe  order  Nj  is  also  equal  to  the  difference 
in  principal  values  of  total  birefringence,  i.e. 


V/hGTS 
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(57) 
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13  Mohr's  circle  of  total  birefringence. 
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When  there  are  no  applied  loads 3  =  ^^2]'/]  “ 

diameter  of  Mohr's  circlej  and  the  value  of  the  fringe  order  seen  in  the 
polariscope,  is  -  N2j^l»  This  is  the  difference  In  the  principal 
values  of  bi ref ringence.  Also  for  this  no-load  cases  ss  v/ill  be  seen^ 
the  isocline  observed  in  the  pel ar1 scope  is  0^0  Thus 3  the  no-load  case 
can  be  used  to  calculate  several  important  parameters  regarding  the 
material's  optical  behavior. 


Redefinition  of  Coefficients 


If  the  following  definitions  are  used 


-3 


f 


2 


f 


12 


N 


R 


(58)  a 

(58)b 

{58)c 

(59) 


eq,  48  becomes 


011  0 
r  i 


2t 

Nj^  cos(29j^)}2  +  s1n(2ep)} 

2  12 


(60) 


With  no  residual  effects,  eq,  60  becomes 


2t 


12 1.2 


^12 


(61) 


This  is  another  form  proposed  by  Sampson  [10].  He  introduced  the  three 
stress-optic  coefficients  f^,  f2,  and  fl2.  These  coefficients  can  be 
determined  for  a  particular  material  by  applying  a  known  stress  and 
observing  the  fringe  order.  Specifically,  if  =  a,  02  =  0,  and 
t^2  =  can  be  computed.  If  =  0,  02  =  a,  and  1-^2  ~  then 

f2  can  be  determined.  If  =  0,  02  =  0,  and  tj2  =  then  fj2  can  be 
calculated.  Alternatively,  knowing  fj  and  f2  and  performing  an  off-axis 
tensile  test,  f^2  <^3n  be  computed. 

For  the  case  of  residual  birefringence,  calibration  of  the  material 
is  more  involved.  With  a  no-load  condition,  and  8^^  can  be  deter¬ 
mined.  Then  with  various  known  stress  states,  f^,  f2,  and  fj2  can  be 
computed. 

For  a  general  case,  knowing  Np,  f^,  f2,  and  the  fringe 

order  can  be  predicted  from  a  given  stress  state,  0^  02,  and  using 
eq.  60.  This  is  exactly  what  will  be  done  later  to  predict  the  fringe 
patterns  in  the  disk. 

In  addition,  for  a  given  stress  state,  the  optical  isocline,  8^, 
can  be  predicted.  Using  eqs.  31e  and  35,  tan(28j)  becomes 


tan(28T) 


2e®  + 

^‘^12M  ^^12M 


T'  /  0  0  'I  ^  R  R  ■v 


IM 


(62) 


With  the  aid  of  eq.  28a-c  this  becomes 


tan(2ey)  = 


2e° 

'^®12M 


+  -  Opp)  sin(28p) 


2RJ 


+  (=1R  -  ®2r)  “s(2eR) 

h 


(63) 


Multiplying  the  numerator  and  denominator  by 
and  47  results  in 


^’'"ave^o 


and  using  eqs.  44 


45 


(64) 


tan(20^)  = 


~  sin(2ej^) 

''’  ■  ^2R^  cos(2ej^) 


Using  eqs.  45  leads  to 


tan(2e^) 


^®66'^12  ^'^IR  “  ^2R^  sin(20j^) 


(a^jj^  -  a2j^)oj^  “  ^^22  ”  ^12^*^2  ^  ^^IR  ~  ^2R^  cos(20j^) 


(65) 


The  definitions  of  eqSo  58  and  59  produce  the  final  form,  namely. 


2 

2  yii  +  N|^  sin(20j^) 

tan(20-r)  = - — - 

T'  0.  Op 

y—  -  -j—  +  Np  cos(20j^) 


(66) 


This  equation  relates  the  optical  isocline  to  the  residual  effects,  the 
material  constants  fj^,  f2,  and  f22s  and  the  stress  state.  The  optical 
isocline  is  the  parameter  observed  in  the  polariscope.  It  is  important 
to  note  that  the  optical  isocline  is  a  function  of  the  stress  level! 
Doubling  each  stress,  for  example,  changes  0j.  The  mechanical  isocline, 
as  0^  is  often  called,  can  be  related  to  stress  from  eq.  66  by  setting 
residual  effects  equal  to  zero.  The  result  is 


tan(20^) 


(67) 


In  the  presence  of  residual  birefringence,  0^  v/ill  never  be  observed. 
Its  existence  can  only  be  hypothesized.  In  the  absence  of  residual  ef¬ 
fects,  the  mechanical  and  optical  isocline  coincide.  The  mechanical 
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isocline  is  not  a  function  of  stress  level.  Doubling  each  stress  has  no 
effect  on  tan(2e^). 

Finally,  with  no  applied  stress  eq.  66  indicates 

Np  sin(20p) 

cos(?6-^T  '  <2®r)  (68) 

or,  as  stated  earlier, 

0j  =  6p.  (69) 

The  observed  isocline  with  no  applied  load  is  the  residual  isocline, 

®R* 

With  the  material  calibrated  for  a  given  stress  state,  the  observed 
isocline  can  be  predicted  using  eq.  66.  This  prediction,  along  with  the 
prediction  of  Nj,  will  be  done  for  the  disk. 

Determining  Stresses  from  Optical  Data 

To  this  point  the  equations  have  been  manipulated  so  that  if  the 
stress  state  is  known,  the  optical  response  of  the  material  can  be  de¬ 
termined.  Specifically,  eq.  60  relates  the  observed  fringe  order  with 
the  stress  state  and  residual  birefringence.  Equation  66  relates  the 
observed  isocline  to  the  stress  state.  To  determine  the  stress  state 
from  the  optical  data,  the  equations  must  be  used  differently.  As  with 
isotropic  photoelasticity,  an  auxiliary  condition  such  as  the  shear  dif¬ 
ference  method  or  a  finite-difference  form  of  the  two  plane-stress 
equilibrium  equations  is  needed  to  completely  determine  the  stress 
state.  However,  photoelastic  data  from  polariscope  observations  pro¬ 
vides  part  of  the  needed  information.  There  are  several  ways  to 
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determine  just  what  information  is  available  from  a  polariscope  in  the 
way  of  the  state  of  stress  in  fiber-reinforced  composites„  The  approach 
taken  here  starts  with  eq.  37. 

In  eq„  37,  subtract  the  2nd  equation  from  the  1st  and  multiply  the 
3rd  equation  by  2.  The  result  is,  with  a  slight  rearrangement 

(qil  -  "  ^^^22  ■  ‘^12^'^2  ""  ■  ®2M^  ‘  '  ®2M^ 

^%6^12  "  ^®12M  "  ^®12M  * 

Using  eqSo  28a-c  and  31a-c  in  the  right  hand  sides  lead  to 


(ill”  ^21^'^!  ~  ^^22”  *^12^*^2  ^^IT”  ®2T^  cos(29|^)  -  ‘'Os(29|^) 

(72) 

"  *^®1T  “  ®2t)  sin(2e^)  -  -  e^,^)  sin(20j^)  (73) 

It  is  obvious  at  this  point  that  multiplying  any  dielectric  tensor 

quantity  by  - ^ —  defines  a  similar  birefringent  tensor  quantity. 

2Xn  e 
ave  0 

Thus  multiplying  both  sides  of  eqs.  72  and  73  by - j - and  using  eqs, 

2xn  e 

46,  47,  and  55-59  results  in  ° 


and 


-  Ny  cos(29-j-)  -  Npj  cos(2e|^) 


^12  ^T 

-y—  =  ^  sin(20y)  -  —  sin(2e|^)  . 


(74) 


(75) 


For  isotropic  photoelasticity,  with  no  residual  effects. 
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=  c  =  =  f^2’ 

and 

Nr  =  0  =  Or.  (77) 

C  being  the  material's  stress-optic  coefficient.  Then  eqs.  74  and  75 
simplify  to 

-  Oy  =  C  N  cos(29)  (78) 

=£jSin(2e)  ,  (79) 

where  the  x-y  system  has  been  substituted  for  the  1-2  material  system. 
These  are  the  classical  equations  of  isotropic  photoelasticity. 

Equations  74  and  75  indicate  that  to  determine  the  stress  state 
from  the  optical  data,  the  difference  in  the  principal  values  of  residu 
al  birefringence,  Nr,  and  the  direction  of  residual  birefringence,  6r, 
must  be  known.  As  stated  previously,  both  these  quantities  can  be 
determined  in  a  polariscope  with  a  no-load  condition. 


Transformation  of  Equations  to  Global  x-y  bystem 

In  reality,  it  is  not  always  convenient  to  reference  the  stresses, 
the  residual  effects,  and  the  material  constants  to  the  material  system. 
Equations  74  and  75  can  be  transformed  to  a  global  x-y  system  which 
makes  an  angle  of  ()>  with  the  1M-2M  system  as  shown  previously.  The 
angle  (}>  is  positive  as  measured  counter-clockwise  from  +x  axis  to  the 
+1M  axis.  With  this  definition  of  <{>,  the  stresses  in  the  material  sys¬ 
tem  are  related  to  the  stresses  in  the  x-y  system  as  follows: 
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COS^(ti 


si 


2sin(j)  cos(j) 

-  2sin4)  cos<j) 
sine})  costfi  sinij)  coscj)  (cos2(})  -  sin^^i) 


sin2(|) 


cos2(i> 


.  (80) 


Substituting  for  and  ^nd  75  results  i 


in 


cos2(j>  a,  +  sin^<b  a  +  2sin(i)  cosi  t 
X  ^  y  xy 

fj 


sin^ij)  a  +  cos^if)  a  -  2sin(j)  cos^i  t 


xy 


and 

-sin(}i  cos(j)  a  +  sintj;  cosij)  a  +  (cos2<j)  -  sin2(j))T 


'Ij  cos(2ej)  (81) 
\  sin(2e^) 


xy 


'12 

—  sin(26j)  -•  —  sin(20j^) 


Gathering  terms  results  in 


(82) 


rCOS^(j)  sin2(}i>|  rsin^ij)  cos^cl)', 

'•~q - fJ~J°x  - Tp^y  ^ 


■1  .  1 


(83) 


2(7^  +  7^)  sin(|)  cosef)  =  N-j.  cos(29^)  -  cos  (20,^)  , 


and 


12 


12 


12 


xy 


"2“  sin(20j)  -  Y~  si0(20|^)  . 


(84) 


If  the  following  definitions  are  made: 
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1  _  rCos^<{)  sin2(}),i 
C  “  ^  f  "  f  '' 

^1  h  ^2 

1  _  i-cos^ip  sin2(j)>| 
c  ~  I  f  “  f  J 

^2  ^2  h 


=  2(-^  +  ■^)  sin({)  cos<t) 

U3 


sin<))  coS(| 
fl2 


1  _  cos^cj)  -  sin^(j) 


then  the  stress-optic  equations  become 


OCT 

Q —  -  “  N-j-  cos(20j)  -  Nj^  cos{20p) 


^  — ~  sin(20y)  -  ^  sin(20|^) 

4  5 


Couple  these  with  the  plane  stress  equilibrium  equations 


(85)a 

{85)b 

(85)c 

(85)d 

{85)e 


ax  *  87  =  ° 

(88)a 

^"^xy  ^  ^®y  _  g 
ax  ay  ’ 

(88)b 

and  a  complete  determination  of  the  stresses  is  possible.  Since  (!>,  fj, 
fp,  and  fi2  are  known,  -  C5  are  known.  Note  that  the  quantities  ob¬ 
served  in  the  polariscope,  Nj,  N;^,  0-|-,  and  0|^,  enter  directly,  indepen¬ 
dent  of  the  x-y  system.  With  certain  situations,  I/C4  and/or  I/C5  might 
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be  zeroo  This  implies  €4  =  00  and  C5  =  if  eqs.  86  and  87  are  being 
implemented  on  a  computer,  it  is  more  convenient  to  define 


C  =  — 
1  C 


1 


r  -  1 

^2  "  iq 


r  -  1 
*=3 


r  “ 
^4 


c  =  — 

^5  C: 


so  that  the  governing  equtions  become 


(89)a 

(89)b 

(89)c 

(89)d 

(89)e 


Cl  ~  ^2  ^3  ■'xy  “  cos(20j)  -  cos(2e,^)  (90) 


«  Nj  Np 

^4  ■  ^4  ‘^y  ^5  "^xy  "  T~  S''ri(2ej)  -  sin(2eji).  (91) 

This  circumvents  having  transformed  stress-optic  properties  becoming 

infinite.  It  is  interesting  to  note  that  when  6  =  45°,  Cr  =  0  and  the 

b 

equation  usually  associated  with  computation  of  the  shear  stress,  eq. 

91,  contains  no  information  about  the  shear.  Also  note  that  for  certain 
materials  the  values  of  fi,  fp,  and  4)  may  combine  to  produce  zero  values 
for  Cl  or  Cg, 
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In  a  later  section,  eqs,  86  and  87  or  90  and  91  will  be  used  to 
determine  the  stresses  in  a  disk  using  the  photoelastic  data  and  a 
finite-difference  representation  of  eqs.  88a  and  b. 


53 


CALIBRATION  OF  THE  MATERIAL 


Uniaxial  Tests 

From  the  previous  section  it  is  clear  that  to  predict  Nj  and  9j  for 
a  given  stress  state ^  the  material  properties  f25  f|2s.  %  and 
must  be  known,  Converselyj  to  predict  the  stress  state  from  Ny  and 
these  same  material  properties  must  be  known.  In  the  present  work^  fol¬ 
lowing  the  approach  of  previous  investigators  [eog,  lllj  an  uniaxial 
loading  of  three  specimens  was  initially  used  to  determine  fj,  f|2. 
Residual  birefringence  was  not  even  considered.  Thus  N[^  and  g^  were  im¬ 
plicitly  set  to  zero  by  using  equations  which  didn't  account  for  residu¬ 
al  effects^,  for  example  eqs,  61  and  67„  Three  calibration  specimens 
were  cut  from  the  same  unidirectional  material  that  was  cut  into  a  disk 
for  later  testing.  The  calibration  specimens  were  cut  as  follows:  one 
parallel  to  the  fibers  {0°)s  one  perpendicular  to  the  fibers  (90°),  and 
one  at  45°  to  the  fiber  direction.  The  material  used  in  these  experi¬ 
ments  was  manufactured  by  the  Illinois  Institute  of  Technology  Research 
Institute  under  contract  with  the  NASA-Langley  Research  Center  [11], 

The  material  had  a  volume  fraction  of  fibers  of  55-60%,  This  was  deter¬ 
mined  by  the  acid  ingestion  method  with  small  pieces  from  the  uniaxial 
specimens.  The  specimens  were  strain  gaged,  to  determine  elastic  prop¬ 
erties,  and  then  subjected  to  uniaxial  loading.  The  fringe  order  was 
observed  as  the  load  level  was  varied.  For  the  tensile  specimens,  at  a 
given  load,  the  entire  specimen  should  have  been  in  the  same  state  of 
birefringence.  Bending  effects  and  other  eccentricities  were  minimized 
and  felt  to  be  negligible.  Yet,  during  the  uniaxial  tests  it  was  not 
possible  to  have  the  entire  specimen  in  what  was  felt  to  be  a  uniform 
state  of  birefringence.  It  was  slightly  lighter  or  darker  toward  the 
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top  of  the  specimen  or  to  one  side,  never  100%  uniform.  This  was  par¬ 
ticularly  true  for  the  higher  load  levels.  Two  observers  were  used  to 

minimize  value- judgement  biases  in  deciding  at  what  load  level  the  spe¬ 

cimen  appeared  to  exhibit,  say,  the  first  fringe.  Observations  were  made 
while  loading,  while  unloading,  on  different  days,  and  for  different 
specimen  orientations  in  the  loading  frame.  At  the  time  the  testing  was 
being  conducted,  it  was  felt  local  variations  in  volume  fraction  or  in¬ 
dex  of  refraction  properties  of  the  fibers  and  matrix  caused  the  slight 
spatial  non-uniformities.  This  idea  was  reinforced  somewhat  by  the  fact 
that  independent  of  orientation  in  the  loading  frame,  the  same  particiu- 
lar  region  of  the  specimen  always  exhibited  slightly  different  birefrin¬ 
gence  effects.  For  the  high  load  levels  in  the  0°  specimen  and  for  all 
load  levels  in  the  90°  and  45°  specimen,  the  material  was  such  that  the 

fibers  were  clearly  at  a  different  level  of  birefringence  than  the 

matrix. 

Figure  14  shows  the  0°  specimen  with  no  load  in  the  lightfield 
polariscope  arrangement.  The  specimen  was  0.625  in.  wide,  about  12  in. 
long,  and  0.090  in.  thick.  In  the  photograph  the  strain  gages  are 
visible,  as  are  the  spots  of  adhesive  securing  the  lead  wires. 

Excessive  strain  gage  adhesive  caused  the  opaqueness  in  the  region 
around  the  strain  gages.  There  was  a  long  hair-sized  flaw  along  the 
length  of  the  specimen  at  the  bottom  end.  This  is  visible  in  the 
photograph.  This  flaw  did  not  appear  to  affect  birefringence  in  the 
region  around  it.  Light  centerlines  scribed  on  the  specimen  are  visible 
between  the  two  strain  gages.  The  intersection  of  these  lines  was 
chosen  as  the  point  at  which  to  count  the  fringe  number.  As  stated 
earlier,  the  birefringent  effect  was  not  spatially  uniform  on  the 
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specimen  and  so  for  consistency  the  fringe  count  was  always  evaluated  at 
this  same  specific  point.  In  the  figure  the  central  part  of  the  speci¬ 
men  appears  lighter  than  the  upper  and  lower  portions.  This  is  because 
the  polari scope  light  source  and  the  specimen  centerline  were  aligned  on 
the  optical  axis  of  the  polariscope  and  the  light  source  was  a  bit  too 
intense.  This  is  not  related  to  the  aforementioned  slight  lack  of  uni¬ 
formity  of  birefringence  in  the  specimen.  Figure  15  shows  the  0®  speci¬ 
men  in  the  light  field  at  what  was  judged  to  be  a  fringe  order  of  1/2. 
Figure  16  shows  the  1st  fringe  in  the  light  field.  Generally,  the  eye 
discerns  spatial  variation  in  light  intensity  differently  than  photo¬ 
graphic  film.  Because  of  this,  except  for  the  aforementioned  difference 
in  intensity  in  the  central  portion  of  the  specimen,  the  birefringence 
of  the  specimens  in  the  photographs  presented  so  far  has  appeared  quite 
uniform  and  is  representative  of  the  uniformity  until  high  load  levels 
were  reached.  Figure  17  shows  fringe  order  4  in  the  dark-field  polari¬ 
scope.  The  lack  of  uniformity  of  the  birefringence  is  evident  in  this 
photograph.  Bundles  of  fibers  could  be  clearly  distinguished  and  there 
was  a  region  in  the  specimen  at  the  lower  left  which  always  had  a  mark¬ 
edly  different  level  of  birefringence.  This  is  the  lighter  region  in 
the  photo.  If  this  had  been  due  to  bending  effects,  there  would  have 
been  an  accompanying  light  region  at  the  top  of  the  test  specimen.  The 
concern  with  lack  of  uniform  birefringence  was  as  follows:  With  a  spa¬ 
tially  varying  stress  state  there  is  a  spatial  variation  of  birefrin¬ 
gence,  If  the  stress  state  is  spatially  uniform  and  the  photoelastic 
sensitivity  of  the  material  varies  spatially,  then  there  is  also  spatial 
variation  in  birefringence.  If  both  the  state  of  stress  and  the  photo¬ 
elastic  sensitivity  vary  spatially  in  an  unknown  manner,  it  would  be 
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difficult,  if  not  impossible,  to  extract  information  regarding  the 
stresses  from  the  fringe  pattern.  Figure  18  shows  the  4th  order  fringe 
in  the  light  field.  The  lighter  region  at  the  lower  left  of  the  speci¬ 
men  in  fig.  17  appears  slightly  darker  in  the  situation  of  fig.  18. 

Figure  19  shows  the  fringe  order  as  a  function  of  stress  for  the  0° 
specimen.  The  experimental  data  are  represented  by  the  solid  circles. 

The  results  of  several  tests  on  different  days  using  two  observers  are 
lumped  together.  If  it  is  assumed  that  there  was  no  initial  birefrin¬ 
gence  in  the  material,  then  when  =  0,  Nj  is  assumed  to  be  zero.  For 
the  case  of  no  initial  birefringence,  the  hypothesis  regarding  the  ef¬ 
fect  of  stress  on  the  dielectric  tensor  ultimately  resulted  in  eq.  61. 
With  02  =  t:j^2  ”  ®  relation  between  N  and  o^  is  a  simple 

straight  line  through  the  origin.  Assuming  no  initial  birefringence, 
i.e,  including  the  point  o^  =  0,  N  =  0,  a  least-square  linear  fit  of  the 
data  from  the  0°  specimen  produced  the  solid  line  in  fig.  19.  This  line 
has  a  slope  of  0.121  x  lO"^  fringe/psi  and  an  intercept  of  0.08  fringe. 
The  slope  led  to  a  value  for  f^  of  744  psi/fringe/in.  If,  as  was  done 
after  seeing  these  data  along  with  the  data  from  the  90  specimen,  it  is 
assumed  that  there  was  some  initial  birefringence,  then  when  Oj  -  0,  Ny 
is  not  necessarily  zero.  This  nonzero  value  of  Ny  at  “  0  is  really 
The  value  of  Nr  could  be  estimated  from  the  data  of  fig.  19  by  not 
using  the  points  =  0,  Ny  =  0  in  the  least-squares  analysis.  Not  in¬ 
cluding  this  data  point  resulted  in  the  linear  least-squares  fit  dotted 
line  in  fig,  19.  The  slope  of  this  line  is  0.119  x  10“^  fringe/psi  and 
it  has  an  intercept  of  0,13  fringe.  The  slope  led  to  fj  =  756  psi/fringe 
/in.  The  calibration  thus  depended  on  whether  or  not  initial  birefrin¬ 
gence  was  assumed.  It  should  be  pointed  out  that  if  at  the  onset  of 
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calibration  residual  birefringence  is  suspected,  the  Tardy  method  can  be 
used  to  determine  the  level. 

Figure  20  shows  the  longitudinal  stress-longitudinal  strain  rela¬ 
tion  for  the  0°  specimen  while  fig,  21  shows  the  lateral  strain-longitu¬ 
dinal  strain  relation  for  that  specimen.  From  least-square  fits. 

Young's  modulus  in  the  fiber  direction,  Ej ,  for  the  specimen  was  5,4  x 
10^  psi  and  the  major  Poisson's  ratio,  0,3, 

Testing  the  90°  specimen  was  more  difficult  than  testing  the  0° 
specimen.  The  ultimate  tensile  strength  of  the  material  in  the  matrix 
direction  was  low  compared  to  the  tensile  strength  in  the  fiber  direc¬ 
tion,  In  addition,  there  was  some  variation  in  the  ultimate  tensile 
strength  in  the  matrix  direction  from  one  point  to  the  next  in  the  ma¬ 
terial,  Thus  it  was  not  always  possible  to  obtain  the  same  maximum 
fringe  count  in  replicate  90°  specimens,  A  90°  specimen  was  tested  in 
tension  until  failure  and  then  a  portion  of  the  failed  specimen  was  tes¬ 
ted  in  compression.  Figure  22  shows  the  set-up  used  for  testing  the  90° 
specimen  in  compression.  The  loading  frame  was  a  hand-operated  screw 
device.  The  1000  lb,  load  cell  used  to  measure  the  compressive  load  is 
shown.  The  specimen  was  strain  gaged  and  it  was  restrained  from 
buckling  by  sandwiching  it  between  two  pieces  of  plexiglas.  The  speci¬ 
men  was  0,625  in,  wide,  3,6  in,  long,  and  0,90  in,  thick.  There  were 
holes  in  the  plexiglas  restraints  so  that  the  restraints  were  not  in 
contact  with  tne  strain  gages.  The  plexiglas  was  not  being  compressed 
even  though  the  specimen  was  and  so  the  fringes  observed  in  the  specimen 
through  the  plexiglas  were  the  true  fringes  in  the  90°  specimen. 

Figure  23  shows  the  no-load,  lightfield  state  of  bi refringence  in 
the  specimen.  Figure  24  shows  the  1/2  order,  lightfield  fringe 
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condition.  Figure  25  illustrates  the  darkfield  state  of  birefringence 
at  a  fringe  order  of  about  1.25.  The  spatial  nonuniformity  of  the  bire¬ 
fringence  is  quite  evident  here,  with  nonuniformity  both  widthwise  and 
lengthwise.  As  one  moves  from  the  top  of  the  specimen  toward  the  bot¬ 
tom,  the  image  of  the  birefringence  is  darker  near  the  centerline  of  the 
specimen  than  it  is  near  the  edges.  This  indicates  the  center  of  the 
specimen  was  at  a  different  level  of  birefringence  than  the  edges. 

About  one-half  way  down,  just  above  the  lower  (lateral)  strain  gage,  the 
birefringence  suddenly  becomes  more  nearly  uniform  with  width.  This 
phenomenon  is  felt  to  be  related  to  inhomogeneities  in  the  material 
rather  than  inhomogeneities  in  the  stress  state.  In  the  90°  direction, 
the  orthotropic  photoelastic  material  crushes  quite  easily.  Thus  any 
deviations  from  flatness  of  the  ends  of  the  specimen,  where  the  compres¬ 
sive  loads  are  transmitted  to  the  specimen,  would  be  eliminated  by  the 
application  of  load.  In  addition,  in  the  90°  direction  end  effects  at¬ 
tenuate  quite  rapidly  due  to  the  low  stiffness.  Thus  any  deviations 
from  flatness  or  other  end  effects  would  be  very  unlikely  to  cause 
abrupt  changes  in  the  widthwise  uniformity  of  birefringence  at  the  cen¬ 
ter  of  the  specimen.  The  halos  aound  the  hexagonal  bolts  were  due  to 
stresses  induced  in  the  plexiglas  when  the  holes  were  drilled.  With  the 
situation  depicted  in  fig.  25  it  is  evident  why  at  times  it  was  diffi¬ 
cult  to  determine  exactly  what  fringe  order  was  being  observed.  For 
this  90°  specimen,  the  fringe  order  at  the  center  of  the  specimen  was 
used  as  the  fringe  order  for  calibration. 

Figure  26  illustrates  the  relation  between  applied  stress  and 
fringe  order  for  the  90°  specimen.  Note  both  tensile  and  compressive 
stresses  are  shown.  As  mentioned  previously,  Nj  is  being  considered  as 


a  positive  quantity  and  so  fig.  26  should  actually  be  plotted  as  shown 
schematically  in  fig.  27.  However,  for  purposes  of  conveniently  using 
least-square  techniques  and  for  strictly  asthetic  reasons,  plots  such  as 
in  fig,  26  will  be  used. 

If  it  is  assumed  that  there  was  no  residual  birefringence,  then 
when  02  =  0>  the  fringe  order  is  zero  (see  eq.  61  with  =  0)« 

The  solid  line  in  fig.  26  shows  the  least-square  calibration  with  the  no 
residual  stress  assumption.  This  line  has  a  slope  of  0.237  x  10 
fringe/psi  and  intercept  of  0.09  fringe.  The  slope  led  to  a  value  of  f2 
=  380  psi/fringe/in.  If  it  is  assumed  that  residual  birefringence 
existed,  the  least-square  data  fit  produced  the  dotted  line.  This  line 
has  a  slope  of  0.227  x  10"^  fringe/psi  and  an  intercept  of  0.15  fringe. 
From  this  slope  f2  =  396  psi/fringe/inch.  Reiterating  an  earlier  dis¬ 
cussion  to  avoid  ambiguity,  the  solid  line  was  obtained  by  using  the 
points  02  =  0,  Nj  =  0  in  the  least-square  analysis.  The  dotted  line  was 
obtained  by  omitting  those  data  in  the  least-squares  analysis.  The 
former  approach  tends  to  force  the  line  through  the  origin  when  in  fact 
it  should  not  go  through  the  origin. 

With  both  figs.  19  and  26  in  hand,  residual  birefringence  was  sus¬ 
pected.  With  fig.  19  alone  it  was  not  clear  whether  the  nonzero  inter¬ 
cept  on  the  N  axis  was  due  to  the  data  scatter  or  due  to  residual  bire¬ 
fringence.  The  nonzero  intercepts  of  both  figures  were,  however,  strong 
arguments  for  residual  effects.  If  there  were  residual  effects,  were 
figs.  19  and  26  consistent?  To  determine  this,  consider  eq.  60. 

With  02  =  becomes 

+  '^R  cos(29,^)}2  +  (N^  sin(29[^)}2  .  (92) 


73 


With  =  0>  eq.  60  becomes 


With  an  actual  polariscope  measurement  it  was  quickly  established  that 
for  the  uniaxial  material,  6j^  was  very  close  to  zero.  Thus  eqs,  92  and 
93  become 


Figures  28a  and  b  illustrate  the  N  vs,  a  relations  for  these  two  equa¬ 
tions.  The  presence  of  residual  birefringence,  with  the  residual  iso¬ 
cline  aligned  with  the  fibers,  results  in  shifts  in  the  calibration  re¬ 
lations.  The  0°  calibration  is  shifted  in  a  direction  opposite  to  the 
shift  of  the  90°  specimen.  The  experimentally  observed  shifts  in  the 
two  specimens  were  consistent  with  the  notion  of  residual  birefringence 
being  present,  particularly  the  sense  of  the  shifts. 

Figures  29  and  30  illustrate  the  elastic  properties  of  the 
90°  specimen.  Figure  29  indicates  that  Young's  modulus  in  the  matrix 
direction  ,  E2,  was  1.78  x  10®  psi  and  fig.  30  indicates  the  minor 
Poisson's  ration  was  0.094. 

A  third  tensile  specimen,  with  the  fibers  at  45°  relative  to  the 
tensile  direction,  was  used  to  complete  the  calibration.  The  data  from 
this  third  specimen  could  be  used  in  several  ways.  The  use  depended  on 
whether  or  not  residual  effects  were  felt  to  be  important.  Figure  31 
shows  the  no-load  state  of  birefringence  in  the  45°  specimen  in  the 
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light  field.  Back-to-back  strain  gage  rosettes  were  used  to  measure  the 
elastic  response  of  the  specimen.  Figure  32  shows  the  lightfield  1/2- 
order  state  of  birefringence  of  the  45°  specimen.  This  was  a  relatively 
low  load  level  and  yet  the  fibers  caused  striations  in  the  birefrin¬ 
gence.  This  could  have  been  due  to  very  high  shear  stresses  at  the 
fiber-matrix  interfaces.  Figure  33  shows  the  45°  specimen  in  the  dark 
field  at  what  was  judged  to  be  the  1st  fringe  order.  As  can  be  seen,  it 
was  quite  difficult  to  determine  what  constituted  a  specific  level  of 
birefringence.  For  this  specimen  the  striations  were  present  at  all 
load  levels.  The  specimen  did,  however,  go  through  lighter  and  darker 
states  of  'striated'  birefringence.  The  lightest  or  darkest  states,  as 
the  load  was  varied,  were  used  to  define  the  fringe  order. 

For  the  45°  specimen  the  combined  stress  state  was  related  to  the 
applied  stress,  o,  by  the  following: 


_  0 
®2  ■  7 


(96)a 


_  a 

^12  ■  7 


(96)b 


If  residual  effects  were  ignored,  these  stresses  could  be  substituted 


into  eq,  61  to  yield 


The  slope  of  the  fringe  order  vs,  applied  stress  relation,  N  vs  a,  can 
be  used  along  with  previously  obtained  values  of  fj^  and  f2  to  determine 
the  value  of  fi2*  34  shows  the  fringe  order-applied  stress  rela¬ 

tion  for  the  45°  specimen  and  also  indicates  a  least-square  straight 


Figo  31  The  45°  tensile  specimen  in  light  field  with  no  load. 
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line.  Since  it  was  assumed  in  fig.  34  that  residual  effects  were  not 
important,  the  points  a  =  0,  Nj  =  0  were  used  as  data  in  the  least- 
square  scheme.  The  slope  of  the  line,  0.288  x  10^  fringe/psi,  resulted 
in  a  value  of  312  psi/fringe/in  for  fi2* 

If  it  is  assumed  residual  effects  were  important  in  the  45°  speci¬ 
men,  eq.  60  must  be  used.  It  is  more  convenient  to  square  eq.  60  to 
obtain 


where  eq. 

96  has  been 

used  along  with  the 

fact  that  =  0. 

If  A,  B, 

and  C  are 

defined  as 

A  ’  It  (f-  -  * 

12 

(99)a 

®  ^  “  T 

K  Tj 

-)  . 

(99)b 

C  =  , 

(99)c 

then  eq.  98  is  of  the  form 


=  Aa2  +  Bo  +  C  .  (lOO)a 

2 

If  the  experimental  measurements  of  Nj  and  o  are  plotted  as  N-j-  vs.  a, 
and  a  least-squares  parabola  of  the  form 

N-j.  =  ao^  +  3o  +  Y  (lOO)b 

is  fit  to  the  data,  then  A,  B,  and  C  can  be  estimated  from 
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A  »  a 


(lOl)a 

(lOl)b 

(lOl)c 


B  «  3 
C  »  Y  , 


If  fj  and  f2  have  been  determined  from  previous  calculations,  then  a  can 

2 

be  used  to  determine  f]^2»  The  value  of  y  should  give  a  value  of  con¬ 
sistent  with  the  value  obtained  in  0°  and  90°  tests.  The  value  of  B 

should  be  consistent  with  the  values  of  N[^,  f]^,  and  f2.  Figure  35  shows 

2 

the  relation  between  Ny  and  a.  The  least-squares  parabola  has  the  fol¬ 
lowing  coefficients, 

a  =  0.75  X  10”^(fringe/psi )2  (102)a 

$  =  0.45  X  10“^(fringe)2/psi  (102)b 

Y  =  0.049  (fringe)^  (102)c 

This  value  of  a  led  to  a  value  for  fj2  of  335  psi/fringe/in.  The  value 
of  Y  predicted  a  value  of  Nj^  =  0.22.  This  value  of  N|^  was  higher  than 
the  value  of  Np  computed  from  the  0°  or  90°  specimens.  However,  there 
was  much  more  difficulty  in  determining  the  level  of  birefringence  in 
the  45°  specimen  and  so,  as  seen  in  figs.  34  and  35,  there  was  much 
scatter  to  the  Nj  -  o  data.  Hence,  stress-optic  properties  obtained 
from  this  specimen  may  not  be  as  reliable  as  properties  taken  from  the 
0°  and  90°  specimens.  Unfortunately,  with  this  approach  to  calibration, 
fl2  must  be  determined  from  this  45°  specimen. 

In  addition  to  the  value  of  Nr  being  different  than  previous  calcu¬ 
lations,  the  value  of  3,  eq.  102b,  was  inconsistent.  The  quantity  3 


85 


approximates  B  in  eq.  100  (see  eq.  99b).  Since  >  f2,  B  will  be  a 
negative  number.  The  value  of  3  from  the  least-squares  approximation 
was  positive.  The  reason  for  this  inconsistency  was  not  evident.  It  is 
not  felt  to  be  due  entirely  to  the  data  scatter  responsible  for  the 
inconsistent  value  of  Np. 

Figure  36  shows  a  plot  of  the  shear  stress  x^2  ^  function  of  the 

sheer  strain  shear  strain  was  determined  from  the  rosette  ga¬ 

ges  while  the  shear  stress  was  computed  using  eq.  96b.  For  some  reason 
the  strain  gage  readings  were  erratic  and  so  the  elasticity  data  had 
scatter.  A  least-square  straight  line  through  the  t  -  y  data  yields  a 
value  of  6^2  570,000  psi.  The  longitudinal  stress-longitudinal 

strain  ralation  for  this  45°  specimen  resulted  in  a  Young's  modulus, 
i.e.  E45,  of  1.56  X  10^  psi. 

Table  1  summarizes  the  optical  properties  and  Table  2  summarizes 
the  elastic  properties  of  the  material  as  obtained  from  the  uniaxial 

tests. 

Four  Point  Bend  Tests 

After  it  was  hypothesized  that  residual  birefringence  was  present 
in  the  material,  it  was  decided  to  use  four-point  bending  tests  to  op¬ 
tically  calibrate  the  material.  There  were  several  reasons  for  this. 
First,  it  was  felt  the  stress  gradients  in  a  bending  specimen  would  re¬ 
sult  in  sharp,  more  well  defined  fringes.  It  would  thus  be  easier  to 
judge  the  presence  of,  say,  the  first  fringe  order.  Using  the  flexure 
formula  and  determining  the  distance  from  the  neutral  bending  axis  to  a 
particular  fringe,  the  stress  levels  producing  that  fringe  order  could 
be  accurately  determined.  Second,  with  residual  birefringence  in  the 


87 


Table  1 


Optical  properties  of  unidirectional  material 
as  determined  from  uniaxial  tests 


Specimen 

Type 


No  Initial 

Birefringence  Assumed 


Initial 

Birefringence  Assumed 


(1) 


0° 

90° 

45° 


fj  =  744  psi/fringe/in 
f2  =  380  psi/fringe/in 
fl2  =  312  psi/fringe/in 


fj  =  756  psi/fringe/in 
Nr  =  0.13 

f2  =  396  psi/fringe/in 
Nr  =  0.15 

fj^2  =  335  psi/fringe/in 
Nr  =  0.22^2^ 


(1)  9r  measured  to  be  zero 

(2)  Nr  and  3  not  consistent  with  previous  computations 
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Table  2 


Elastic  properties  of  unidirectional  material 
as  determined  from  uniaxial  tests 


Specimen  Type 


Elastic  Properties 


0° 

90° 

45° 


=  5.4  X  10^  psi 
^2  =  1*78  X  10^  psi 
^12  “  0.57  X  10®  psi 


V 22  “  0.3 
V22  =  0.09 


E45  =  1.5  X  10®  psi 


90 


material,  the  fringes  would  not  be  symmetrically  disposed  on  either  side 
of  the  bending  neutral  axis.  The  'optical  neutral  axis'  would  be  off¬ 
set  from  the  conventional  neutral  axis,  i.e.  the  centerline  of  the  bend¬ 
ing  specimen.  This  off-set  could  be  used  to  determine  the  residual 
fringe  order  N|^.  Finally,  one  test  could  be  used  for  both  tensile  and 
compressive  stresses.  With  a  known  applied  bending  moment  and  using 
Tardy  compensation  to  determine  fractional  fringe  orders,  considerable 
stress-optic  data  could  be  obtained  from  the  four-point  bending  arrange¬ 
ment. 

All  the  previously  discussed  uniaxial  specimens  were  used  as  four- 
point  bending  specimens.  These  specimens  were  reused  as  bending  speci¬ 
mens  for  two  reasons.  First,  it  conserved  valuable  photoelastic  materi¬ 
al.  Second,  it  insured  that  variability  of  optical  properties  from  lo¬ 
cation  to  location  on  the  panel  would  not  contribute  to  any  differences 
in  the  calibration  constants  between  the  uniaxial  tests  and  the  four- 
point  bend  tests.  There  was  one  disadvantage  in  using  the  uniaxial 
specimen.  These  particular  specimens,  as  stated  before,  were  approxim¬ 
ately  0.625  in.  wide.  Wider  specimens  would  have  resulted  in  more  ac¬ 
curacy  in  resolving  distances  from  the  neutral  bending  axis.  Actually, 
photographs  of  the  four-point  bend  test  fringe  patterns  were  enlarged  to 
aid  in  determining  distances  from  the  neutral  axis.  Thus  the  size  of 
the  specimen  was  not  important  from  that  standpoint.  However,  using  a 
wider  specimen  as  the  basis  for  scaling  would  have  resulted  in  more  res¬ 
olution  than  using  a  narrower  one. 

Figures  37-39  show  the  three  four-point  bending  tests.  Figure  37 
shows  the  0°  specimen.  In  this  case  the  fibers  were  parallel  to  the 
neutral  bending  axis  (horizontal  in  the  figure).  Figures  38  and  39, 
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Fig.  38  Four-point  bend  test  of  90°  specimen. 


Fig.  39  Four-point  bend  test  of  45°  specimen 


respectively,  show  the  90°  and  45°  specimens.  With  the  90°  orientation, 
the  fibers  were  perpendicular  to  the  neutral  axis.  The  90°  specimen  was 
somewhat  shorter  than  the  other  two  since  it  was  the  remains  of  a  failed 
tensile  specimen.  Recall,  the  90°  tensile  specimens  were  tested  in  ten¬ 
sion  to  failure  and  then  one  of  the  remaining  pieces  was  tested  in  com¬ 
pression  with  the  buckling  restraints. 

With  the  specimens  as  shown  in  figs.  37-39,  the  portion  of  the  beam 
above  the  neutral  axis  was  in  compression.  For  the  0°  specimen, 


(103) 


where  z  was  measured  positive  downward  from  the  neutral  axis  and  M  and  I 
have  their  usual  meaning.  Equation  94  provides  the  relation  between 
and  Nj,  i.e. 


Using  various  values  of  moment,  looking  at  several  z  locations,  and 
using  Tardy  compensation,  sufficient  data  could  be  generated  to  reliably 
determine  f^  and  N[^.  Figure  40  shows  the  relation  between  and  Nj. 

The  maximum  stress  levels  used  in  the  four-point  tests  were  not  neces¬ 
sarily  identical  to  the  stress  levels  in  the  uniaxial  tests  because  in¬ 
terest  here  was  in  residual  birefringence.  The  slope  of  the  line  in 
fig.  40  is  0.106  x  10"^  fringe/psi  and  the  intercept,  N[^,  is  0.11.  The 
slope  led  to  a  value  of  fj  of  849  psi/fringe/  inch.  This  is  11%  higher 
than  the  value  of  the  f^  using  the  tensile  tests,  assuming  an  initial 
birefringence. 
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For  the  90°  specimen. 


°2 


Mz 

I 


(104) 


and  eq„  95  provides  the  relationship  between  and  N-j-,  i.e. 


Figure  41  shows  the  relation  between  Nj  and  02*  The  slope  of  the  rela¬ 
tionship  is  0e234  X  10“^  fringe/psi  and  the  intercept  is  -0.07  fringe. 
The  slope  resulted  in  a  value  of  f2  of  385  psi/fringe/inch.  This  is 
about  3%  lower  than  the  value  obtained  in  uniaxial  tests. 

For  the  45°  specimen,  eqs.  96  and  98  are  applicable,  where  in  those 
equations 


a  = 


Mz 

I  * 


(105) 


Figure  42  shovifs  the  relation  between  Nj  and  a.  A  least-squares  parabola 
of  the  form  eq.  100b  was  fit  through  the  data.  The  coefficients  were 


a  -  0.837  X  10“^  (fringes)2/(psi)2  (106)a 

g  =  0,34  X  10“^  (fr1nges)2/psi  (106)b 

y  =  0.0451  fringes  .  (106)c 

Using  the  values  of  fi  and  f2  from  the  0°  and  90°  bending  tests,  the 
value  of  a  led  to  a  value  of  fj^2  319  psi/fringe/in.  This  is  about 
10%  lower  than  the  value  of  fj^2  obtained  from  the  uniaxial  tests.  The 
value  of  y  led  to  a  residual  fringe  value  of  0.21.  The  sign  of  3  was 
again  inconsistent  with  the  other  data.  Table  3  summarizes  the  optical 
calibration  from  the  four-point  bending  tests. 
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Fig.  42  Stress-optic  relation  for  45°  specimen  from  four-point  bend  test. 


iable  3 

Optical  properties  of  unidirectional  niaterlal 
as  determined  from  four-point  bend  tests. 


=  849  psi/fringe/in 
Nr,  -  Ooll 

f2  =  385  psi/fringe/in 
Nr  =  0o07 

=  319  psi/fringe/in 
Nr  =  0,21^2) 


9~  measured  to  be  zero 

&  not  consistent  with  previous  coaiputaticns 


0° 

90“^ 

45° 
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In  summary,  it  is  felt  that  the  four-point  bending  tests  yielded  a 
more  accurate  calibration.  The  ambiguity  concerning  fringe  definition 
in  the  uniaxial  tests  did  not  occur  in  the  bending  tests.  Furthermore, 
having  tensile  stresses  and  compressive  stresses  at  the  same  time  was 
felt  to  be  valuable.  Since  the  transition  between  tension  and  compres¬ 
sion  is  associated  with  determining  the  value  of  N|^,  having  a  single 
testing  fixture  and  simultaneously  determining  the  value  of  Nj  associa¬ 
ted  with  tension  and  the  value  associated  with  compression  should  reduce 
the  variance  in  the  determination  of  N|^. 
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FORWARD  TESTING:  GLOBAL  BEHAVIOR 


With  the  material  calibrated,  eqs.  60  and  66  could  be  used  to  pre¬ 
dict  the  isochromatic  and  isoclinic  fringes  in  a  disk,  assuming  the 
stresses  were  known.  Alternatively,  using  eqs.  74  and  75,  in  addition 
to  a  third  condition,  the  stresses  could  be  computed,  assuming  the  iso¬ 
chromatic  and  isoclinic  fringe  patterns  were  known.  The  former  has  been 
referred  to  here  as  forward  testing  while  the  latter  has  been  referred 
to  as  backward  testing.  The  forward  testing  is  discussed  in  this  and  in 
the  next  section 

The  comparing  of  the  predicted  fringe  patterns  with  the  experiment¬ 
ally  observed  fringe  patterns  was  done  using  two  approaches.  Both  com¬ 
parisons  were  done  using  an  orthotropic  disk  in  diametral  compression. 
The  first  comparison  was  done  by  looking  at  the  disk  as-a-whole.  This 
is  discussed  in  this  section  and  is  referred  to  as  global  behavior.  The 
second  approach  is  described  in  the  next  section  and  is  referred  to  as 
local  behavior.  To  study  the  global  behavior,  photographs  of  fringe 
patterns  in  the  disk  were  compared  with  computer-generated  fringe  pat¬ 
terns.  These  computer-generated  fringe  patterns  were  based  on  eq.  60 
and  the  elasticity  solution  of  Appendix  B.  This  solution  is  based  on 
the  complex  variable  approach  of  Lekhnitskii  [12].  The  stress-optic 
data  required  in  eq.  60  was  taken  from  the  tensile  and  four-point  bend¬ 
ing  tests  and  the  stresses,  0^,  ©2,  and  tj2»  were  determined  from  the 
elasticity  solution.  Figure  43  shows  the  set-up  used  to  load  the  disks 
in  diametral  compression.  The  disk  used  was  3,12  in.  in  diameter  and 
0.090  in,  in  thickness.  Shown  in  the  photograph  is  the  plexiglas  sand¬ 
wiching  plates  used  to  prevent  the  disk  from  buckling  under  the 
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compressive  load.  This  loading  fixture  is  the  fixture  described  in 

conjunction  with  the  discussion  of  the  90°  compression  calibration 

specimens o 

One  dilemma  in  the  comparison  of  theory  and  experiment  was  the 
choice  of  the  numerical  values  of  the  stress-optic  parameters^  particu¬ 
larly  The  calibration  tests  produced  several  values  of  ranging 
from  0o07  to  0,22.  Realistically,  the  0,22  value  from  Table  1  and  the 
0,21  value  from  Table  3  were  unreliable.  As  previously  indicated,  they 
were  not  consistent  with  the  other  calibration  data  from  the  same  speci¬ 
men,  It  was  not  clear  a  simple  average  of  the  various  values  would  be  a 

valid  choice  for  the  value  of  N[^,  A  value  of  0,15  was  chosen,  somewhat 

arbitrarily,  for  This  certainly  was  a  good  approximation  to  the 

values  measured  from  the  tests.  The  values  of  the  other  stress-optic 
coefficients  were  as  follows: 

f;^  =  849  psi/fringe/in, 

f2  =  385  psi/fringe/in, 
fl2  =  319  psi/fringe  in. 

These  data  are  from  Table  3, 

Figure  44  shows  a  typical  computer-generated  fringe  pattern  for  the 
disk.  For  the  case  shown,  the  fringe  order  is  Nj  =  0,5  and  the  fibers 
are  30°  relative  to  the  loaded  diameter.  Figures  45,  48,  51,  54,  and  57 
show  the  comparisons  between  the  observed  fringe  patterns  in  the  disk 
and  the  computer-generated  fringe  patterns.  The  figures  show  the  fringe 
patterns  for  the  disk  with  the  load  at  various  angles  relative  to  the 
fibers.  In  all  figures  the  primary  fringe  order  is  Nj  =  0,5,  Other 


(107)a 

(107)b 

(107)c 
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Fig.  44  Computer-generated  f 
compression,  Nj=0,5 


fringes  are  visible  near  the  application  of  the  load„  These  fringes 
were  not  studied  because  of  the  lack  of  ideal  loading  and  lack  of  accur¬ 
acy  in  the  elasticity  solution  in  this  region„  Visible  in  the  photo¬ 
graphs  of  the  disk  is  crazing  damage  on  the  disk  edge  appearing  as  dark 
spots  on  the  edge,,  When  the  load  was  applied;,  the  disk  material  in  the 
immediate  vicinity  of  the  load  application  point  would  craze  and  thus 
become  opaque  to  light  transmission.  At  first,  the  load  was  introduced 
into  the  disk  with  a  rounded  contact  point.  The  radius  of  curvature  of 
the  rounded  contact  was  0,125  in.  This  small  radius  was  used  to  simu¬ 
late  point  loads.  However,  the  damage  induced  in  the  disk  was  severe. 
Thus  the  flat  contacts  seen  in  the  figure  were  used.  Still  some  crazing 
v/ould  always  occur,  especially  when  loading  the  disk  perpendicular  to 
the  fibers. 

Figure  45  shows  the  observed  and  predicted  fringe  pattern  when  the 
load  was  applied  along  the  fibers.  The  load  level  was  500  lb.  Certain¬ 
ly  the  characteristics,  such  as  the  vertical  bow-tie  shape,  correlated 
well.  On  the  one  hand,  it  is  quite  remarkable  how  good  the  correlation 
'was.  The  state  of  stress  in  the  disk  was  quite  complex.  The  stress- 
optic  law  is  a  hypothesis  and  the  calibration  of  the  stress-optic  coef¬ 
ficients  was  done  with  simple  tests.  The  three  main  ingredients  seemed 
to  combine  well  to  lend  credibility  to  the  whole  notion  of  orthotropic 
photoelasticity.  On  the  other  hand,  a  closer  scrutiny  of  the  results 
shows  the  predicted  horizontal  distance  from  the  center  of  the  disk  to 
the  Nj  =  0,5  fringe  was  less  than  observed.  Scaled  to  the  disk’s  actual 
size,  the  distance  from  the  center  to  the  fringe  was  predicted  to  be 
0,33  in.  From  the  photograph  the  distance  was  closer  to  0,40  in.  In¬ 
terestingly  enough,  if  it  was  assumed  the  residual  fringe  value,  N^,  was 
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ted  and  observed  fringe  pattern 
Nt=0.5  (predicted=white  dots) 


45  Comparison  between  predic 
fibers  aligned  with  load 


0,07 5  one-half  the  value  being  assumed ,  the  distance  prediction  was  verv’ 
close  to  the  value  observed.  Figure  46  shows  the  predicted  fringe  pat¬ 
tern  as  a  function  of  three  levels  of  residual  bi  ref  ri  ngence ,  =  0-  i'b 

=  0.07,  and  Nf^  =  0,15,  The  level  of  residual  birefringence  obvicusly 
makes  a  difference  in  the  fringe  pattern  generated  in  the  disk,  Figir  r- 
47  shows  another  interesting  aspect  of  the  effects  of  residual  bire¬ 
fringence,  The  figure  shows  the  predicted  effects  of  simply  re'.'ersl nc; 
the  load  and  subjecting  the  disk  to  diametral  tensiun.  ’-or  no  resiurai 
birefringence,  the  sign  of  the  load  has  no  influence  on  tlie  observed 
fringe  pattern.  With  residual  bi refringance,  the  sign  of  the  load 
influences  the  fringe  pattern  produced. 

Figure  48  shows  a  comparison  between  predicted  and  obser-vad  fringe 
pattern  in  the  disk  when  the  load  was  oriented  30°  relative  to  the  fi 
bers.  The  load  level  was  325  lb,  (With  the  sign  convention  adapted  in 
Appendix  B,  this  corresponds  to  a  =  -30°,)  The  fibers  caused  the 
streaking  in  the  photo  and  also  contributed  to  the  lack  of  sharp  fri,ige 
definition  in  some  regions  of  the  disk.  This  obscuring  was  the  same  as 
observed  in  the  testing  of  the  uniaxial  specimens  described  earlier. 

The  distance  from  the  disk  center  to  the  fringe  location,  along  the  di¬ 
ameter  perpendicular  to  the  fibers,  was  slightly  cverpredi cted  by  the 
theory.  Figure  49  shows  the  effect  of  various  levels  of  residual  bire¬ 
fringence  on  the  predicted  fringe  patterns.  The  differences  are  subtle 
but  there  are  slight  differences.  When  shear  stresses,  in  the  in  -  am 
coordinate  system,  dominate  the  stresses,  small  changes  in  residual  bi- 
ref  ringent  effects  become  less  important.  This  is  because  for  this  liid- 
terial  Qr  =  0  and  only  the  term  involving  a-|  and  a.^  ir,  eq,,,  60  is  affect¬ 
ed  by  varying  Mr,  Figure  50  shows  the  effect  of  a  load  reversal  on  the 
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predicted  fringe  pattern  for  this  30°  orientation.  Again,  it  is  clear 
that  simply  reversing  the  load  direction  has  an  impact  on  the  fringe 
pattern . 

Figure  51  shows  the  predicted  and  observed  fringe  patterns  for  the 
case  of  the  fibers  making  a  45°  angle  with  the  diametral  compressive 
load.  The  load  level  was  325  lb.  (This  orientation  corresponds  to  a  = 
+45°.)  Again  the  fibers  tended  to  spread  the  fringe  but  the  pattern  was 
quite  close  to  the  predicted  configuration.  A  close  scrutiny  of  the 
distance  from  the  disk's  center  to  the  half-order  fringe,  along  a  diame¬ 
ter  perpendicular  to  the  fibers,  showed  the  predicted  value  to  be  quite 
close  to  the  observed  value.  Figures  52  and  53  show  the  effects  of 
residual  level  and  load  reversal,  respectively. 

Figures  54-56  shows  similar  comparisons  for  the  case  with  the 
fibers  60°  relative  to  the  load  direction.  Again  the  load  level  was  325 
lb.  (This  corresponds  to  a  =  +60  in  the  nomenclature  of  Appendix  B). 

The  spreading  of  the  fringe  by  the  fibers  was  quite  severe.  However, 
the  computer-generated  fringe  pattern  shows  good  correlation  between  the 
dots  of  the  predicted  fringe  pattern  and  the  center  of  the  fringe. 

Figure  57  shows  the  observed  and  predicted  fringe  pattern  for  the 
situation  when  the  load  was  90°  relative  to  the  fibers.  Here  the  load 
level  was  520  lb.  The  correlation  for  this  case  was  the  poorest  of  all 
cases.  Using  a  residual  birefringence  of  0.07  made  the  correlation 
better.  Figure  58  shows  that  the  predicted  fringe  pattern  for  this  ori¬ 
entation  of  load  and  fibers  is  fairly  sensitive  to  the  level  of  residual 
birefringence.  As  expected,  Fig.  59  shows  that  with  residual  bire¬ 
fringence  the  sign  of  the  load  is  important  in  determining  the  nature  of 
the  observed  fringe  pattern. 
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Effect  of  residual  birefringence  on  the  load  reversal  fringe 
pattern,  fibers  45°  relative  to  load,  Nt=0,5, 


birefringence  on  predicted  fringe 
to  load,  Nj=0.5. 


Fig.  57  Comparison  between  predicted  and  observed  fringe  pattern, 

fibers  perpendicular  to  load,  Nj=0.5  (predicted=white  dots). 
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g.  59  Effect  of  residual  birefringence  on  the  load  reversal  fringe 
pattern,  fibers  perpendicular  to  load,  Nj=0,50, 
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In  conclusion^  it  appears  that  the  idea  of  orthotropic  photoelasti- 
crcy  is,  on-the-v/hole,  a  valid  concepto  However,  fringe  definition  can 
definitely  be  a  problem.  In  the  cases  here,  the  'answer'  was  known  and 
so  it  is  easy  to  conclude  the  fringes  w'ere,  more  or  less,  where  they 
should  have  been.  However,  for  the  60°  orientation  (fig,  54),  at  the 
two  o'clock  and  eight  o'clock  locations  on  the  disk  (high  noon  being  at 
the  upper  load  application  point),  it  may  not  have  been  possible  to 
sketch  in  the  fringes  if  it  were  not  known  a  priori  where  they  were  sup¬ 
posed  to  be  located.  This  is  a  serious  concern  in  the  practical  appli¬ 
cation  of  orthotropic  photoelasticity.  This  prompted  the  series  of 
tests  described  in  the  following  section,  namely,  a  closer  look  at  the 
predicted  and  observed  fringes. 


FORWARD  TESTING:  LOCAL  BEHAVIOR 


Whereas  forward  testing  of  the  global  behavior  indicated  an  overall 
comparison  between  theory  and  experiment,  a  more  detailed  comparison  was 
of  value.  For  this  reason  another  series  of  comparisons  between  theo¬ 
retical  and  experimental  fringe  data  was  made.  Specifically,  the  total 
isochromatic  fringe  order,  Nj,  and  the  optical  isocline,  0^,  were  com¬ 
pared  point-by-point  across  a  diameter  perpendicular  to  the  load  and 
across  a  diameter  perpendicular  to  the  fibers.  These  comparisons  were 
made  for  diametral  loads  at  various  angles  relative  to  the  fibers. 

These  comparisons  are  shown  in  figs.  60-75.  The  figures  show  the  exper¬ 
imentally  measured  values  and  their  theoretical  predictions.  On  each 
figure  there  are  several  theoretical  predictions,  each  prediction  corre¬ 
sponding  to  a  different  level  of  assumed  residual  birefringence,  Np. 

The  three  levels  are  the  same  as  in  the  previous  section,  =  0,  0.07, 
and  0.15.  The  isochromatic  fringe  predictions  were  based  on  eq.  60,  the 
stress-optic  data  for  the  material,  and  the  complex  variable  elasticity 
solution.  The  optical  isocline  predictions  were  based  on  eq.  66,  the 
stress-optic  data  for  the  material,  and  the  elasticity  solution.  On 
each  figure  the  horizontal  axis  shows  the  nondimensional  distance  across 
a  diameter  as  measured  from  the  left  edge  of  the  specimen.  Measurements 
were  only  made  from  one  edge  of  the  disk  to  the  center.  The  vertical 
axes  show  the  isochromatic  or  isoclinic  measure,  depending  on  the  fig¬ 
ure.  On  the  figures,  the  open  squares  represent  the  experimentally  mea¬ 
sured  values.  The  solid  circles  represent  the  theoretical  predictions 
for  the  case  of  no  residual  birefringence  {N[^  =  0).  The  solid  triangles 
represent  the  predictions  for  the  case  of  N^j  =  0.07  and  the  solid 
squares  represent  the  theoretical  values  for  the  case  =  0.15. 
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Figure  60  shows  the  comparison  of  isochromatic  fringe  order  for  the 
case  when  the  loading  was  aligned  with  the  fibers.  The  load  level  used 
was  407  Ibo  Obviously  the  diameter  perpendicular  to  the  load  and  the 
diameter  perpendicular  to  the  fibers  were  the  same  for  this  case.  As 
can  be  seen,  starting  at  the  left  edge  of  the  disk  and  moving  toward  the 
center,  the  experimentally  measured  isochromatic  fringe  count  decreased 
to  practically  zero  and  then  increased  to  a  value  slightly  greater  than 
one-half  at  the  center.  Independent  of  the  assumed  level  of  initial  bi¬ 
refringence,  the  theoretically  predicted  fringe  count  also  approached 
zero  somewhere  along  the  radius.  The  initial  level  of  birefringence 
only  affected  where  it  went  to  zero.  For  no  initial  birefringence  (sol¬ 
id  circles)  the  fringe  count  dropped  to  zero  at  x/D  =  0.15.  Experimen¬ 
tal  measurements  (open  squares)  showed  minimial  fringe  value  at  x/D  = 
0.21.  A  residual  level  of  birefringence  of  0,07  caused  the  predicted 
fringe  level  to  drop  to  zero  at  x/D  =  0.23.  The  zero  fringe  was  pre¬ 
dicted  to  occur  at  x/D  =  0.3  v/hen  the  residual  birefringence  was  0.15, 
Except  for  near  the  left  edge,  the  experimental  observation  correlated 
well  with  the  predictions  based  on  low  (Nj^  between  0  and  0,07),  as  op¬ 
posed  to  high  (N[^  =  0,15),  levels  of  residual  birefringence. 

Figure  61  shows  the  predicted  and  observed  optical  isocline  for  the 
case  of  the  load  aligned  with  the  fibers.  This  figure  is  quite  inter¬ 
esting.  In  reality  all  symbols  should  be  plotted  at  6  =  0°  or  6  =  -90°. 
To  illustrate  an  important  point,  the  symbols  have  been  shifted  slightly 
off  the  0°  and  -90°  values.  It  was  the  rapid  jump  of  the  isocline  from 
0°  to  -90°  as  one  moved  in  along  a  radius  that  was  important.  The  x/D 
location  at  which  the  optical  isocline  jumped  from  0°  to  -90°  was  the 
same  location  at  which  the  total  isochromatic  fringe  count  dropped 
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fringe  order. 


x/D 


Fig.  60  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  diameter  perpendicular  to  load, 
fibers  aligned  with  load. 
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isoclinic  angle. 


Fig„  61  Comparison  between  predicted  and  observed  isoclinic 
angle  along  diameter  perpendicular  to  load,  fibers 
aligned  with  load. 


to  zero.  This  location  can  be  thought  of  as  an  optical  isotropic  point. 
At  such  a  point  the  total  isochromatic  fringe  order  is  zero  and  the 
isoclinic  is  indeterminant.  Referring  to  eq.  60,  the  second  squared 
bracketed  term  is  zero  along  the  diameter.  There  is  no  shear  stress, 
and  the  residual  isocline,  0„,  is  zero.  With  a  proper  combination 

IC  r< 

of  normal  stresses,  and  a^,  the  first  squared  bracketed  term  can  be 
made  zero.  This  was  precisely  the  case  at  the  various  x/D  locations. 

As  can  be  seen  from  eq.  60,  the  value  of  Np  affects  this  term  going  to 
zero.  Examining  eq.  66,  if  the  first  and  second  terms  of  eq.  60  are 
zero,  eq.  66  becomes 

tan  29j  =  .  (108) 

This  leads  to  i ndeterminent  values  for  6^.  It  should  be  noted,  not  only 
does  the  value  of  Mr  affect  the  location  of  the  isotropic  point,  but  the 
load  level  does  also.  Since  Nr  is  fixed  and  and  can  be  varied  by 
varying  the  load  level,  the  physical  location  of  the  optical  isotropic 
point  can  be  controlled.  This  was  not  explored  in  this  study. 

From  the  above  discussion,  it  seems  that  it  may  be  possible  to  de¬ 
termine  Nr  by  comparing  the  theoretical  and  experimental  isotropic  point 
location  along  the  diameter.  Since  the  isotropic  point  in  the  experi¬ 
ment  was  located  between  the  predicted  isotropic  point  for  no  initial 
birefringence  and  the  case  Nr  =  0.07,  the  assumed  lower  levels  of  resid¬ 
ual  birefringence  seem  to  correlate  better  with  experiment  than  does  the 
assumed  higher  level. 

Figures  62  and  63  show  further  isochromatic  and  isoclinic  data  for 
the  disk  with  the  load  aligned  with  the  fibers.  Since  the  diameter 
perpendicular  to  the  load  and  the  diameter  perpendicular  to  the  fiber 
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Fig.  62  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  30°  diameter,  fibers  aligned  with 
load. 
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Fig.  63  Comparison  between  predicted  and  observed  isocline 
angle  along  30°  diameter,  fibers  aligned  with  load. 


were  tha  sa:n83  data  was  studied  along  another  clanetar.  This  diameter 
y/as  arbitrarily  picked  to  be  a  diameter  vshlch  made  a  30'’  angle  with  the 
h  G  r 1 2  0  n  t  a 1  d1 ame  taro 

Moving  fronn  the  edge  of  the  disk  tov/ard  the  center  along  the  30° 
dlameoer  produced  trends  in  the  data  similar  lo  the  horizontal  diameter 
oota,  The  isochromatic  fringe  number  decreased  and  than  Increased  to 
about  0o5.  The  minimum  value  of  the  traoreti cal ly  predicted  Isochromat- 
ic  varieds  depending  on  the  assumed  value  of  Initial  bi  ref^'i  ngenceo  For 
no  assLimed  initial  b1  refri rigence  (solid  circles).,  the  1  sochromati c 
fringe  value  dropped  to  zero  at  x/D  =  0„i2,  The  experimental ly  measured 
1 sochromati c  never  dropped  to  zero.  The  same  was  true  for  the  theoreti¬ 
cal  predictions  based  on  assumed  initial  birefringence  values  of  0.07 
and  0.15. 

Figure  63  shows  the  theoretically  predicted  and  the  experimentally 
maasiired  optical  isocline  pararnater  for  the  30°  diameter.  There  'were 
vast  ciifferenceSj  depending  on  the  assumed  level  of  initial  birefring¬ 
ence,  between  the  predicted  and  the  measured  values.  The  experimentally 
measured  isocline  varied  smoothly  from  about  5°  at  the  outer  edgCj 
through  0°  at  x/0  =  0.06 5  to  -85°  at  the  canter  of  the  disk.  The  theor¬ 
etical  cases  of  Ni^  =  0,07  and  Np^  =  0,15  also  varied  smoothly  but  the 
case  Np  =  0  exhibited  a  sharp  jump  at  x/D  =  0.13,  This  was  the  same  lo¬ 
cation  that  Nj  =  0  for  this  theoretical  case.  This  behavio''  was  similar 
to  the  jump  behavior  in  the  optical  Isocline  along  the  horizontal  diame¬ 
ter  for  tir's  fiber/load  orientation.  Agairia  except  for  the  outer  edge. 


theoreti ca ! 


predictions  for  the  case  Mp  =  0,07  seemed  to  correlate 


the  best  with  experimental  observations.  The  1 sochromati c  predictions 
almost  coincided  and  the  isocline  predictions  were  close. 
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Recall,  there  was  concern  over  choosing  the  value  of  Np  from  the 
tensile  and  compression  test  data.  There  were  values,  though  unrelia¬ 
ble,  of  Np  as  high  as  0.21.  There  was  strong  evidence  from  the  tensile 
and  compression  data  that  N[^  was  in  the  range  0.07  -  0.15.  The  value  of 
0.15  was  chosen  from  the  test  data  of  the  0°,  90°,  and  45°  specimens. 

Yet  the  experimentally  observed  global  and  local  behavior  correlated 
better  with  theoretical  predictions  obtained  by  using  a  lower  value  of 
assumed  residual  birefringence,  specifically  N[^  =  0.0  or  0.07.  This 
trend  toward  good  correlation  with  0  <  Np  <  0.07  in  the  global  and  local 
tests  contradicted  the  findings  of  the  0°,  90°,  45°  tests.  The  contra¬ 
diction  is  still  a  puzzle. 

Figures  64  and  65  show  the  isochromatic  and  the  isoclinic  behavior 
along  a  horizontal  diameter  for  the  case  of  the  fibers  rotated  30°  from 
the  applied  diametral  load.  The  observed  isochromatic  and  isoclinic  be¬ 
havior  was  similar  to  the  other  two  cases  discussed.  The  experimentally 
observed  isocline  changed  from  negative  values  to  positive  values  at  x/D 
of  about  0.08.  The  theoretical  prediction  for  the  case  Nf^  =  0  showed  a 
very  sharp  jump  in  isocline  at  this  x/D  location.  At  the  same  time,  the 
isochromatic  prediction  for  this  level  of  birefringence  practically 
dropped  to  zero  at  this  same  location.  Throughout  the  testing  of  all 
cases  presented  in  the  figures  in  this  section,  when  the  isocline  exper¬ 
ienced  a  rapid  change  with  spatial  location,  the  isochromatic  fringe 
number  was  going  to  a  minimum  value.  This  was  true  with  the  experiment¬ 
al  observations  and  all  theoretical  predictions.  It  seemed  that  the 
closer  the  isochromatic  fringe  approached  zero,  the  more  dramatic  was 
the  jump  in  the  isoclinic. 
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Fig.  64  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  diameter  perpendicular  to  load, 
fibers  30°  relative  to  load. 
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Fig.  65  Comparison  between  predicted  and  observed  isoclimc 
angle  along  diameter  perpendicular  to  load,  fibers 
30®  relative  to  load. 
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Several  additional  conments  should  be  made  regarding  the  data  in 
figs.  64  and  65.  In  the  range  0.2  <  x/D  <  0.5,  the  assumed  level  of  re¬ 
sidual  birefringence  did  not  affect  the  predicted  i sochromati c.  All 
predictions  coincided.  Because  of  this,  in  that  same  range  along  the 
diameter,  the  correlation  between  predicted  and  observed  isochromatics 
didn't  seem  to  favor  any  particular  level  of  Np.  However,  the  isoclinic 
figure  showed  reasonable  correlation  with  the  N|^  =  0.07  case,  particu¬ 
larly  the  characteristics  near  the  change  in  sign  in  the  isoclinic. 

Figures  66  and  67  show  more  isochromatic  and  isocline  data  for  the 
case  with  the  fibers  oriented  30°  relative  to  the  load.  These  figures 
are  for  the  diameter  perpendicular  to  the  fiber  direction.  Observations 
and  comments  similar  to  the  discussion  of  previous  figures  can  be  made 
regarding  these  two  figures.  As  seen  in  fig.  67,  for  example,  the  ex¬ 
perimentally  observed  isoclinic  experienced  a  sharp  change  at 
x/D  “  0.14.  It  was  at  this  location  Ny  was  a  minimum.  At  first,  the 
sharp  changes  in  isoclines  were  viewed  as  an  undesirable  feature.  Even¬ 
tually,  however,  these  sharp  changes  were  a  welcome  characteristic. 

They  served  as  sort  of  a  datum  for  minimum  isochromatic  fringe  counts 
and  a  point  where  there  was  most  likely  a  change  in  sign  in  the  isoclin¬ 
ic.  Polari scopes,  by  their  nature,  respond  similarly  to  isoclinics  of 
complementary  values,  i.e.,  isoclinics  that  differ  in  value  by  90°, 

When  the  sign  of  the  isoclinic  changes,  or  the  isoclinic  otherwise  goes 
to  zero,  interpretation  of  the  data  can  be  confusing.  The  existence  of 
the  sharp  changes  helped  clarify  the  interpretation. 

Figures  68-75  show  the  isocline  and  isochromatic  data  for  cases 
when  the  fibers  were  at  45°  and  90  to  the  load  direction.  These  cases 
will  not  be  discussed  in  detail.  In  light  of  previous  discussions,  some 


136 


fringe  order. 


0.8 


0.7 


0.6 


0.5 


0.4 


0.3 


0.2 


®  theory,  =  0 

^  theory,  =  0.07 

@  theory, 

□  experiment 


Nn  =  0.15 


0.1 


n 


n 


0.1 


0.2 


□ 


n 


□ 


n 


n 


n 


n 


X  n 


disk  diameter  =  D 


0.3 


I  i 

0.4  0.5 

x/D 


Fig.  66  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  diameter  perpendicular  to  fibers, 
fibers  30°  relative  to  load. 
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Fig,  67  Comparison  between  predicted  and  observed  i sod  ini ( 
angle  along  diameter  perpendicular  to  fibers,  fibei 
30°  relative  to  load. 
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Fig.  68  Comparison  between  predicted  and  observed  isochromatic  fringe 
order  along  diameter  perpendicular  to  load,  fibers  45“  relative 
to  load. 
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F1g.  69  Comparison  between  predicted  and  observed  isoc 
diameter  perpendicular  to  applied  loads,  fiber 
load. 
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Fig.  70  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  diameter  perpendicular  to  fibers, 
fibers  45°  relative  to  load. 
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Fig.  71  Comparison  between  predicted  and  observed  isoclinic 
angle  along  diameter  perpendicular  to  fibers,  fibers 
45°  relative  to  load. 
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Fig.  72  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  diameter  parallel  to  fibers,  fibers 
perpendicular  to  load. 
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Fig.  74  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  30°  diameter,  fibers  perpendicular 
to  load. 
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of  the  observations  and  comments  are  obvious.  The  load  level  used  for 
the  45°  and  90°  orientations  was  350  lb.  It  can  be  seen  from  the 
figures  that  the  experiments  correlated  better  with  the  lower  levels  of 
assumed  residual  birefringence  than  they  did  with  the  higher  levels. 

It  was  of  interest  to  put  these  experiments  of  local  behavior  into 
context  by  repeating  the  experiments  with  a  standard  isotropic  disk.  A 
great  effort  is  made  in  the  manufacture  of  isotropic  photoelastic  mater¬ 
ial  to  eliminate  residual  birefringence.  Thus  the  isotropic  tests  could 
not  be  used  to  check  the  effects  of  residual  birefringence.  However, 
the  tests  could  certainly  be  used  to  check  the  spatial  distribution  of 
isochromatic  and  isoclinic  data.  If  it  is  assumed  the  theoretically 
predicted  data  were  correct,  then  the  experiments  would  evaluate  how 
well  photoelastic  data  was  being  measured.  If  it  was  assumed  that  the 
photoelastic  measurements  were  correct,  then  the  comparison  between  ex¬ 
periment  and  theory  could  be  viewed  as  an  evaluation  of  the  tneory  of 
isotropic  photoelasticity.  The  theory  of  isotropic  photoelasticity  has 
been  verified  time  and  time  again  and  so  experiments  with  the  isotropic 
disk  should  be  looked  upon  as  an  indication  of  the  accuracy  to  which 
isochromatic  and  isoclinic  fringe  data  could  be  recorded  with  the 
pol ari scope. 

Figures  76-79  show  the  comparison  between  the  measured  and  the  pre¬ 
dicted  isochromatic  and  isoclinic  fringe  information  for  the  isotropic 
disk.  The  predictions  are  based  on  the  exact  elasticity  solution  for  an 
isotropic  disk  as  presented  in  Frocht  [13].  The  load  level  used  was  135 
lb.  The  isotropic  disk  material  was  PSM-1.  Figure  76  shows  the 
isochromatic  fringe  numbers  across  a  diameter  perpendicular  to  the  load. 
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Fig.  76  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  diameter  perpendicular  to  load, 
isotropic  disk. 


Figure  77  shows  the  isoclinic  data  for  that  diameter.  The  correlation 
was  very  good.  Figures  78  and  79  show  similar  information  for  a  diame¬ 
ter  oriented  at  30°  to  the  horizontal  diameter.  Again  the  correlation 
was  excellent.  From  these  four  figures,  there  was  no  question  as  to  the 
ability  to  make  accurate  measurements  with  the  polariscope. 

From  these  studies  with  the  orthotropic  and  isotropic  disks,  one 
conclusion  is  certain.  Orthotropic  photoelasticty,  as  it  is  hypothe¬ 
sised  here,  is  not  as  accurate  as  isotropic  photoelasticity.  Although 
the  conclusion  is  straightforward,  it  is  not  so  clear  why  there  is  a 
difference  in  accuracy.  It  can,  however,  be  traced  to  tv/o  sources. 
Possibly,  the  hypothesis  regarding  the  stress-optic  behavior,  given  by 
eq.  36,  is  not  entirely  valid.  The  averaging  of  the  effects  of  the  fi¬ 
bers  and  matrix,  in  a  through-the-thi ckness  sense,  may  not  be  as  accur¬ 
ate  for  optical  behavior  as  it  is  for,  say,  elastic  behavior.  In  addi¬ 
tion,  there  may  be  an  interaction  between  shear  and  elongation  effects. 
That  is,  maybe  eq.  36  should  not  be  orthotropic  in  any  coordinate  sys¬ 
tem.  These  issues  are  directed  at  the  theory  itself.  Perhaps  the  inac¬ 
curacy  comes  in  the  practicing  of  the  theory.  Fringe  definition  in  or¬ 
thotropic  photoelastic  material  is  not  as  good  as  it  is  in  isotropic  ma¬ 
terials.  This  can  be  seen  in  a  number  of  previous  figures,  for  example, 
fig.  54.  Thus  the  discrepancies  between  theory  and  experiment  may  result 
from  not  being  able  to  accurately  locate  the  center  of  a  fringe.  Fringe 
count  was  kept  low  in  the  orthotropic  material,  Ny  -  0.5,  whereas  it  was 
made  high  in  the  isotropic  material,  N  “  5.0.  Due  to  streaking  by  the 
fibers,  the  higher  order  fringes  (greater  than  3  or  4)  are  difficult  to 
define  in  the  orthotropic  material.  In  an  isotropic  material,  higher 
order  fringes  are  narrower  and  more  distinct  than  the  lower  order 
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Fig.  78  Comparison  between  predicted  and  observed  isochromatic 
fringe  order  along  30°  diameter,  isotropic  disk. 
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Fig.  79  Comparison  between  predicted  and  observed  isoclinic 
fringe  order  along  30°  diameter,  isotropic  disk. 
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fringes.  Thus  the  two  materials  respond  in  an  opposite  manner  to  high 
fringe  numbers.  Where  the  stresses  are  high  and  changing  rapidly,  the 
isotropic  material  responds  in  an  ideal  manner.  The  orthotropic  materi¬ 
al  responds  exactly  opposite  to  the  ideal  manner.  This  lack  of  spatial 
resolution  of  fringes  in  the  orthotropic  material  is  felt  to  be  the  ma¬ 
jor  reason  for  the  discrepancies  between  theory  and  experiment.  Also, 
perhaps  a  lack  of  spatially  uniform  residual  birefringence  contributes 
to  the  behavior. 

The  real  issue,  though,  is  how  accurately  can  stresses  be  predicted 
using  orthotropic  photoelasticity.  This  is  the  subject  of  the  next 
section. 
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BACKWARD  TESTING:  PREDICTION  OF  STRESSES 


Historically,  photoelasti ci ty  has  been  used  to  predict  stresses 
from  known  isochromatic  and  isoclinic  fringe  information.  For  a  general 
state  of  stress,  at  each  point  in  a  body  the  two  normal  stresses,  o 

X 

and  a  ,  and  the  shear  stress,  x  ,  are  to  be  determined.  At  each  point 
in  the  body,  the  isochromatic  fringe  count  gives  one  quantity  of  experi¬ 
mental  information  and  the  isoclinic  gives  a  second  quantity  of  experi¬ 
mental  information.  A  third  quantity  or  condition  is  needed  if  three 
stresses  are  to  be  uniquely  determined  at  the  point.  In  isotropic  ma¬ 
terials,  interferometric  measurements  have  been  used  to  provide  informa¬ 


tion  regarding  the  sum  of  principal  stresses.  This  is  a  third  quantity 

and  thus  0,0,  and  x  can  be  determined,  Interferofietri c  measure- 
X  y 

ments  are  not  that  easy  to  obtain  so  experimenters  tend  to  rely  on  other 
relations  among  the  stresses  as  a  third  condition.  One  other  condition 
is  that  the  unknown  stresses  must  satisfy  the  compatibility  equations. 
Compatabi 1 i ty  is  a  second  order  partial  differential  equation.  For  iso¬ 
tropic  materials  the  compatibility  equation  reduces  to  Laplace's  equa¬ 
tion  for  the  quantity  [a  +  a  ].  Since  the  solution  to  Laplace's  equa- 

X  y 

tion  on  simple  domains  such  as  a  circle  or  rectangle  is  well  estab¬ 
lished,  investigators  have  used  the  compatibility  equation  as  the  auxil¬ 
iary  condition.  For  orthotropic  materials,  the  compatibility  equation 
in  terms  of  stresses  does  not  reduce  to  Laplace's  equation. 

Another  condition  often  used  by  photoel asti ci ans  is  that  the  un¬ 
known  stresses  must  satisfy  equilibrium.  Equilibrium  is  represented  by 
two  first-order  partial  differential  equations.  These  two  equations  in¬ 
volve  the  three  stresses  and  by  themselves  cannot  be  solved,  as 
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Laplace's  equation  can,  on  any  domain.  In  addition,  the  two  equilibrium 
equations  plus  the  isochromatic  and  the  isoclinic  fringe  data  result  in 
four  conditions  from  which  to  solve  for  the  three  stresses.  In  this 
situation  the  stresses  are  overdetermined.  To  have  a  unique  solution, 
one  of  the  equilibrium  equations  can  be  simply  discarded.  This  is  done 
in  the  so-called  shear  difference  approach.  If  both  equilibrium  equa¬ 
tions  are  retained  the  stresses  computed  are  the  three  stresses  which 
satisfy  the  four  conditions  in  a  least-squares  sense.  It  is  also  possi¬ 
ble  to  use  the  equilibrium  equations  in  such  a  way  that  in  the  region  of 
interest  in  a  photoelastic  model  there  are  exactly  as  many  unknown 
stresses  as  there  are  equations.  This  approach  was  used  in  this  present 
study  and  it  relies  on  applying  the  equilibrium  equations  in  an  average 

sense. 

Independent  of  how  the  equilibrium  equations  are  used  as  an  auxili¬ 
ary  condition,  the  equilibrium  equations  are  usually  represented  in  fi¬ 
nite-difference  form.  The  partial  derivations  are  approximated  by  using 
central,  forward,  and  backward  differences  in  the  stresses  at  the  gird, 
or  mesh,  points  of  a  finite-difference  grid.  The  partial  differential 
equations  thus  become  algebraic  equations  relating  the  three  stresses  at 
the  mesh  points.  In  isotropic  photoelasticity,  eqs.  78  and  79  are  then 
used  as  two  more  algebraic  equations  at  each  mesh  point.  Since  the  iso¬ 
chromatic,  N,  and  the  isoclinic,  9,  are  known  at  each  point,  eqs.  78  and 
79  represent  algebraic  relations  the  normal  and  shear  stresses  must  sat¬ 
isfy  at  each  point.  The  two  finite-difference  forms  of  the  equilibrium 
equations  and  eqs.  78  and  79  represent  the  four  equations  which  are  used 
to  compute  the  three  stress  components.  A  previous  photoelastic  study 
[14]  of  isotropic  double-lapped,  double-pin  connectors  used  this  four- 
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equation  scheme  to  estimate  stresses,,  This  scheme  was  adapted  to 
orthotropic  photoelastic  materials  by  using  eqs.  86  and  87  in  place  of 
8qs„  78  and  79.  This  is  possible  because  the  form  of  the  equilibrium 
equations,  unlike  the  form  of  the  compatibil ity  equation,  does  not 
depend  on  the  inaterial  properties.  The  finite-difference  representation 
or  the  plane-stress  equilibrium  equation  remains  unaltered.  Appendix  C 
describes  the  scheme. 

To  determine  the  degree  to  which  orthotropic  photoelasticity  could 
be  used  to  predict  stresses,  a  final  series  of  experiments  was  conducted 
with  the  disk  using  the  finite-difference  formulation.  The  disk  was 
loaded  in  diametral  compression  and  the  isochromatic  and  the  isoclinic 
fringe  data  determined  along  various  diameters  of  the  disk.  These  data 
were  used  in  the  finite-difference  formulation  to  calculate  the  three 
stresses  along  these  diameters.  These  stresses  were  then  compared  with 
the  values  computed  using  the  elasticity  solution  of  Appendix  A.  Since 
the  level  of  residual  birefringence  was  an  issue,  the  experimentally 
predicted  stresses  were  computed  using  various  values  of  assumed  residu¬ 
al  birefringence.  The  values  of  assumed  were:  0,  0.07,  and  0.15. 

The  value  of  was  assumed  to  be  0.  For  given  values  of  Nj  and  9^,  the 
experimentally  computed  values  of  a^,  a^,  and  would  be  dependent  on 
Npo  Changing  Nj^  and  changes  the  right-hand  side  of  eqs,  86  and  87. 
Physically,  changes  in  the  right  hand  side  of  the  equation  'adjust'  the 
icinge  data  so  only  those  data  due  to  the  applied  load  generate  numeri¬ 
cal  values  for  the  stresses.  Obviously  if  Nj  =  Np  and  0  =  e  the 
right  hand  side  'adjusts'  to  zero  and  the  system  of  equations  predicts 
that  there  are  no  stresses  due  to  applied  loads. 
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Befors  using  dctual  photoGlastic  data  in  the  finite-difference 

schene,  the  accuracy  of  the  finite-difference  formulation  was  assessed. 

This  was  done  by  using  the  exact  elasticity  solution  and  the  stress- 

optic  law  to  generate  Nj  and  0j  at  each  point  of  the  finite-difference 

mesh.  These  values  of  Ny  and  6^  were  used  in  the  finite-difference 

scheme  as  if  they  were  experimental  data.  The  values  of  a^,  and 

T  were  then  computed.  These  values  of  stress  were  compared  with  the 

xy 

values  of  stress  at  each  mesh  point  as  predicted  by  the  elasticity  solu¬ 
tion  used  to  generate  Nj  and  ®y»  This  exercise  tested  the  accuracy  of 
the  finite-difference  formulation.  In  essence,  if  the  exact  values  of 
Ny  and  8^  were  used,  could  the  finite-difference  scheme  compute  the  ex¬ 
act  value  of  the  stresses?  Generally  the  answer  to  the  question  was  a 
resounding  yes.  For  a  few  cases,  as  will  be  seen,  there  were  some  dif¬ 
ferences  between  the  exact  values  of  and  from  the  elasticity 

solution  and  the  values  from  the  finite-difference  scheme. 

Figure  80  shows  the  finite-difference  mesh  used.  In  the  following 
figures  the  stresses  along  the  middle  line  of  the  three-line  mesh,  i.e. 
the  diametral  line,  are  reported.  The  mesh  consisted  of  0.092  in. by 
0.092  in.  squares.  With  the  disk  3.12  in.  in  diameter,  there  were  18 
grid  points  along  the  disk  radius.  Studies  were  not  conducted  to  deter¬ 
mine  the  effects  of  mesh  refinement  on  predicted  stresses.  Because  of 
the  aforementioned  good  agreement  between  the  elasticity  solution  stres¬ 
ses  and  the  finite-difference  predictions  using  the  exact  values  of  Ny 

and  6y,  such  refinement  studies  were  felt  to  be  superfluous. 

Figures  81-83  show  the  stresses  along  the  diameter  perpendicular  to 
the  load  for  the  case  of  the  fibers  aligned  with  the  load.  The  stresses 
have  been  normalized  by  the  absolute  value  of  the  mean  stress  across  the 


157 


disk  diameter,  !P|/Dt,  where  P  is  the  load,  D  is  the  disk  diameter,  and 
t  is  the  disk  thickness.  (In  the  elasticity  solution  P  is  negative  if 
compressing  the  disk.)  The  normalized  stresses  are  plotted  on  the  ver¬ 
tical  axes  and  the  distance  along  the  diameter  is  plotted  on  the  hori¬ 
zontal  axis.  The  diametral  distance  is  measured  from  the  left  edge  of 
the  disk.  Experimental  measurements  of  Ny  and  9^  were  obtained  for  only 
one-half  the  disk.  The  inset  in  each  figure  indicates  the  fiber  orien¬ 
tation  relative  to  the  load,  and  the  diameter  being  considered. 

Figure  81  shows  variation  of  the  major  compressive  stress,  o^,  with 
diametral  distance.  There  are  several  sets  of  data  on  this  and  on  sub¬ 
sequent  figures.  The  solid  squares  represent  the  predictions  of  the 
elasticity  solution.  The  open  diamonds  represent  the  predictions  of  the 
finite-difference  scheme  when  the  elasticity  solution  values  of  Nj  and 
0j  were  used  as  'experimental'  data.  This  is  the  aforementioned  finite 
difference  check.  In  fig.  81  the  diamonds  are  almost  entirely  obscured 
by  the  solid  squares,  indicating  the  aforementioned  good  accuracy  of  the 
finite-different  approach.  The  slight  differences  that  did  occur  were 
at  the  outer  edge  of  the  disk  where  the  stresses  were  quite  low.  The 
open  squares  represent  the  stress  predictions  when  the  values  of  Nj  and 
Oj,  as  measured  from  the  actual  disk,  are  used  in  the  finite-difference 
scheme  and  the  residual  birefringence  was  assumed  to  be  =  0.15.  The 
open  triangles  indicate  the  experimentally  predicted  stresses  when 
was  assumed  to  be  0.07.  The  open  circles  indicate  the  case  when  was 
assumed  to  be  zero.  As  can  be  seen,  and  as  could  be  expected,  the  level 
of  assumed  residual  birefringence  affected  the  stress  predictions.  The 
close  agreement  between  the  elasticity  solution  (solid  squares)  and  the 
Np  =  0.07  case  was  in  keeping  with  the  findings  of  the  forward  testing 
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Fig.  81  Comparison  of  a  stresses  along  diameter  perpendicular 
to  load,  f ibers'^al  igned  v;ith  load. 


experiments  described  in  the  previous  section.  There,  recall,  with  Nf^  = 
0.07  or  lower,  the  level  of  residual  birefringence  resulted  in  better 
correlation  with  experimental  findings  than  other  levels  of  assumed 
residual  birefringence. 

Figure  82  shows  the  matrix  direction  stress,  a^,  for  this  same 
case.  In  this  figure  the  open  diamonds  are  not  obscured  by  the  exact 
solution's  solid  squares.  It  is  important  to  note  these  stresses  are  an 
order  of  magnitude  less  than  the  stresses  of  the  previous  figure  and 
the  scale  for  the  stress  coordinate  has  been  expanded  to  fill  the  fig¬ 
ure.  Plotting  this  figure  with  the  same  scale  as  used  in  fig.  81  would 
significantly  'squeeze'  the  a  data  together.  What  was  interesting 
about  the  experimental  a  computation  is  that  it  did  not  depend  on  the 

A 

assumed  level  of  residual  birefringence.  This  is  evidenced  by  the  fact 
that  the  open  circles,  squares,  and  triangles  are  coincident.  This 
always  occurred  when  the  x-y  system  coincided  with  the  principal  materi¬ 
al  axes. 

Figure  83  shows  the  shear  stress  calculations,  .  for  the  case 
being  discussed.  Due  to  symmetry  in  geometry,  loading,  and  material 
properties  the  shear  stress  was  zero  across  the  horizontal  diameter. 

The  shear  stresses  shown  in  fig.  83  were  the  values  for  the  upper  mesh 
points  of  the  three-line  finite-difference  grid  of  fig.  80.  These  shear 
stress  were  quite  small  and  really  almost  meaningless.  The  calculations 
could  be  easily  influenced  by  any  inaccuracy  in  measurement  or  computa¬ 
tion.  However,  the  trend  and  magnitude  of  the  theoretically  predicted 
shear  stresss  were  well  matched  by  the  experimentally  determined  values. 
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Figures  84-86  shov'i  the  three  stress  components  along  a  diameter 
ted  30'’  relative  to  the  horizontal.  The  stresses  are  computed  in  a 
inate  system  which  has  its  x  axis  along  this  diameter  and  the  y 
ticn  perpendicul ar  to  this.  The  orientation  of  the  x-y  system  is 


in  the  figures, 

Flqur’e  84  shows  the  variation  of  a.,,  along  this  diameter.  The 

^  O' 

1  designations  ares  will  bes  consistent  v/ith  previous  notation, 
the  stress  calculations  based  on  the  experimental  values  of  Nj  and 
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'r-~'  -r-  r:  V'l 

i  1 


■Cl  li  rs 


earn  to  align  with  any  particular  level  of  assumed  residual 
e,  Hov/ever,  the  experiifiental  calculations  were  more  closely 
the  theoretical  predictions  when  lov;  levels  of  residual 
e  were  assumed  in  the  finite-difference  calculations, 

85  illustrates  the  stress  computations.  Again  the  scale 


rences  bet'ween  this  figure  and  fig,  84  should  be  noted. 

Figure  86  illustrates  the  stress  calculations  For  the  30° 

xy 

tar.  These  stresses  were  on  the  true  diameter  rather  than  for  an 


et  line  as  in  fig,  83,  For  the  majority  of  the  diameter,  the  ex- 
ental  calculations  using  Mp,  =  0,07  coincided  with  the  theoretical 
icity  solution. 

Figures  87--i04  show  the  stresses  for  other  load  orientations  rela¬ 
te  the  fibers.  These  situations  will  not  be  discussed  in  detail 
t  for  cases  w^here  there  were  some  unusual  features.  In  some  cases 
orrelation  between  the  exact  solution  and  calculations  based  on  ex- 
antal  data  was  similar  to  the  reasonable  correlation  for  the  cases 
discu'ssad  in  detail.  In  the  other  case,  correlation  vias  not  as 
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Fig.  87  Comparison  of  a  stresses  along  diameter  perpendicular  to  load, 
fibers  30°  relative  to  load. 
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Fig.  97  Comparison  of  a  stresses  along  diameter  perpendicular  to  fibers, 
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Fig.  101  Comparison  of  t  stresses  just  of  diameter  perpendicular  to  load, 
fibers  perpendicular  to  load. 
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Fig.  102  Comparison  of  a  stresses  along  30°  diameter,  fibers  perpendic¬ 
ular  to  load.  ^ 
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Fig.  103  Comparison  of  stresses  along  30°  diameter,  fibers  perpendic¬ 
ular  to  load. 
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Fig.  104  Comparison  of  x  stresses  along  30°  diameter,  fibers  perpendicular 
to  load. 
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The  shear  stresses  In  fig.  89  showed  a  significant  sensitivity  to 
the  assumed  value  of  M|^.  It  appeared  that  if  a  value  of  =  0.02  or 
0.03  had  been  used,  the  correlation  between  theory  and  experiment  for 
the  shear  stresses  for  that  particular  load  and  diameter  orientation 
would  have  been  good.  The  stresses  for  this  case,  fig.  87,  would 
have  also  correlated  well.  To  keep  the  scatter  of  fig.  89  in  context, 
however,  the  scale  of  fig.  89  should  be  compared  with  the  scale  of  fig. 
87. 

Fig.  90  shows  sensitivity  of  the  major  compressive  stress  component 
to  the  assumed  value  of  Here  again  a  value  of  =  0.02  -  0.03 
would  have  given  excellent  correlation. 

Figure  91  is  the  first  case  to  show  differences  between  the  exact 
elasticity  solution  and  the  finite-difference  calculations  (solid  squares 
and  open  diamonds,  respectively)  when  the  exact  values  of  Nj  and  0.^.  were 
used.  Due  to  the  magnitude  of  the  stress  in  Fig.  91  relative  to  the 
major  stress,  fig,  90,  there  was  no  concern  for  the  differences. 

The  case  shown  in  figs.  93-95  deserves  mention.  Referring  to  eq„ 
85e,  the  value  of  <}>  for  figs.  93-95  is  ^  =  45°.  This  means  I/C5  is  0, 
or  Cg  ^  =0.  Also,  from  eqs.  85a  and  85b,  =  -C2.  This  is  the  special 

case  discussed  in  eqs.  89a-91.  It  required  a  different  computer  pro¬ 
gram,  The  exact  solution  and  the  finite-difference  calculations,  solid 
squares  and  open  diamonds,  were  in  good  agreement  for  this  case.  This 
gave  confidence  in  this  special-case  computer  program.  This  particular 
load/fiber  orientation  resulted  in  poor  correlations  between  theory  and 
experiment,  e.g.  fig.  93. 

The  cases  where  the  load  was  perpendicular  to  the  fibers,  figs.  99- 
104,  offered  some  unusual  comparisons.  The  exact  elasticity  solution 
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and  the  finite-difference  calculations  using  the  exact  values  for  Nj 
and  O-r  deviated  more  in  the  o  stress  than  they  did  for  any  other  case 
studied.  Strangely  enough,  the  and  stresses,  figs.  99,  101,  102, 
and  104  showed  good  comparison  between  the  two  'exact'  cases. 

As  in  the  previous  section  on  forward  testing,  it  was  of  interest 
to  determine  how  closely  the  stresses  could  be  computed  experimentally 
in  an  isotropic  disk.  Figures  105-107  show  the  three  components  of 
stress,  in  an  isotropic  disk,  across  a  diameter  perpendicular  to  the 
load.  The  shear  stress  was  computed  for  the  top  line  of  the  three-line 
finite-difference  grid.  There  was  no  residual  birefringence  in  the 
material . 

Figure  105  shows  the  major  compressive  stress,  o^.  Here  as  in  the 
previous  discussions,  the  solid  squares  are  the  exact  elasticity  solu¬ 
tion  while  the  open  circles  are  the  results  computed  from  the  experimen¬ 
tal  values  of  N  and  9  and  the  finite-difference  scheme.  The  comparison 
was  good.  Exact  values  of  N  and  6  from  the  elasticity  solution  were 
used  to  test  the  finite-difference  calculations  and  the  results  of  that 
calculation  are  displayed  as  open  diamonds.  It  is  impossible  to  see 
them  in  fig.  105  since  they  coincide  with  the  solid  squares. 

Figure  106  shows  the  stress  calculations  while  fig.  107  shows 
the  shear  stress  calculations.  The  comparisons  between  the  exact  solu¬ 
tion  and  the  experimental  calculations  were  good  for  these  stresses.  In 
fact,  they  were  much  better  than  the  comparisons  for  the  minor  stresses 

fa  and  x  )  in  the  orthotropic  material.  Figures  108-110  show  the 
'■X  xy 

stress  comparisons  for  a  diameter  which  was  30°  relative  to  the  diameter 
perpendicular  to  the  load.  The  correlation  was  again  good. 
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Fig.  106  Comparison  of  o  stresses  along  diameter  perpendicular 
to  load,  isotropic  disk. 
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From  the  results  for  the  orthotropic  and  isotropic  disks,  it  is 
clear,  as  it  was  in  the  section  on  forward  testing,  orthotropic 
photoelasticity  does  not  appear  to  be  as  accurate  as  isotropic  photo- 
elasticityo  One  is  faced  with  the  question.  Is  it  the  hypotheses  and 
assumptions  that  are  lacking,  or  is  it  the  practicing  of  the  theory? 

The  lack  of  good  fringe  definition  could  certainly  lead  to  scatter  in 
the  stresses  predicted  from  the  experimental  data.  Uncertainty  in,  or 
simple  miscalculation  of,  the  residual  birefringence  certainly  could 
also  produce  the  poor  results.  As  has  been  seen,  simply  adjusting  the 
level  of  residual  birefringence  can  lead  to  excellent  correlation  be¬ 
tween  the  theoretical  elasticity  prediction  of  stresses  and  the  value 
determined  from  the  experimental  photoelastic  data.  All  evidence  seemed 
to  indicate  the  level  of  residual  birefringence,  and  perhaps  a  lack  of 
spatial  uniformity,  were  responsible  for  the  lack  of  correlation  that 
can  be  achieved  with  an  isotropic  disk.  To  look  at  the  lack  of  correla¬ 
tion  from  one  other  point  of  view,  a  numerical  study  of  the  effects  of 
inaccuracies  in  the  isocline  measurements  was  conducted.  This  study, 
and  the  motivation  for  it,  are  discussed  in  the  next  section. 


EFFECTS  OF  ISOCLINE  ERRORS 
ON  STRESS  PREDICTIONS 

From  the  figures  and  discussion  of  the  fringe  patterns  in  the  or¬ 
thotropic  disk,  e.g.  fig.  51,  it  is  clear  that  fringe  definition  was  a 
concern  in  the  experiments.  The  figures  of  fringe  patterns  shown  previ¬ 
ously  are  illustrations  of  the  isochromatic  fringe.  The  isoclinic 
fringe  have  not  been  illustrated.  Several  isoclinic  fringes,  along  with 
their  computer-generated  predictions,  could  have  been  shown  in  the  sec¬ 
tion  on  forward  testing:  global  behavior.  Since  for  stress  predictions 
local  behavior  was  more  of  a  concern  than  global  behavior,  that  exercise 
was  not  pursued.  However,  the  importance  of  the  isoclines  should  not  be 
under-emphasised.  In  fact,  in  isotropic  photoelasticity,  the  stress 
calculations  are  generally  much  more  sensitive  to  inaccuracies  in  the 
isoclinics  than  to  inaccuracies  in  the  isochromatics.  In  isotropic  pho¬ 
toelasticity,  the  isoclinic  fringes  are  much  harder  to  discern  than  the 
isochromatic  fringes.  The  isoclinics  are  obscured  by  the  isochromatic 
fringes  and  they  are  not  as  sharp.  In  addition,  in  regions  of  low 
stress  gradients  it  is  very  difficult  to  see  the  isoclines.  So  it  was 
with  these  experiment  with  orthotropic  materials.  The  optical  isocline 
was,  at  times,  much  more  difficult  to  see  than  the  isochromatics.  Dis¬ 
persion  by  the  fibers,  regions  over  which  the  isoclines  did  not  change 
rapidly,  and  lack  of  vivid  contrast  between  light  and  dark  regions  of 
the  model,  all  contributed  to  uncertainties  in  the  isoclinic  data.  For 
this  reason,  numerical  studies  were  conducted  to  determine  the  effect  of 

uncertainties  in  the  isocline  data  on  the  predictions  of  a  ,  a  ,  and 

X  y 

T  for  the  various  fiber  and  diameter  orientations.  There  were  strong 

xy 
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indications  that  some  of  the  discrepancies  observed  between  the  exact 
stress  calculations  and  the  computations  based  on  experimental  data  were 
due  to  uncertainties  in  the  residual  birefringence  level.  However^  for 


completeness,  the  effects  of  uncertainties  in  the  isocline  measurements 
were  studied. 

The  Isocline  uncertainty  studies  ware  of  tivo  types.  Both  were 
based  on  the  isoclines  predicted  by  the  exact  elasticity  solution  and 
tl’se  stresS“Optic  laws.  These  were  the  same  data  used  to  generate  the 
open  diamonds  in  figs,  8L-104,  In  the  first  uncertainty  study,  the  the¬ 
oretical  value  of  the  isocline  of  each  mesh  point  in  the  finite-differ¬ 
ence  scheme  vms  increased  by  5°,  This  adjusted  isocline  was  then  used 
v/ith  the  exact  Nj  as  'experimental  data'  in  the  finite-difference  calcu¬ 
lations,  All  stresses  of  all  eight  fiber/diameter  orientations  were 
computed  using  the  skev/ed  isocline.  These  stresses  were  compared  with 
the  stresses  from  the  finite-difference  cal culati ons  using  the  correct 
values  of  8^  and,  of  course,  Nj,  The  results  of  this  study  were  used  to 
determine  hov/  a  5°  error  in  all  isoclines  affected  the  predicted  stres¬ 
ses,  The  value  of  5°  was  chosen  as  a  value  which  could  actually  occur 
if  several  experimental  effects  v;ere  additive.  The  polariscope  being 


used  v/as  a  split  bench  type  and  the  polarizer  and  analyser  were  linked 
together  by  shafts,  universals,  and  other  coupling  mechanisms.  In  addi¬ 
tion,  these  optical  elements  were  rotated  by  means  of  gears.  It  is  en¬ 
tirely  possible  that  backlash  in  the  gears  and  coupling  elements  could 
produce  a  constant  bias  in  the  measurements.  In  addition  to  the  mechan¬ 
ical  bias  in  the  polariscope,  the  alignment  of  the  disks  in  the  polari¬ 
scope  was  another  source  for  a  constant  bias.  It  is  possible  that  the 
fiber  direction  of  the  disk  was  not  aligned  exactly  as  assumed  in  the 
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polariscope.  There  is  a  0°  direction  in  the  polariscope  but  there  were 
no  cross-hairs  or  other  alignment  aids  to  insure  that  the  direction  of 
polarization  of  the  optical  elements  were  aligned  at  a  particular  angle 
relative  to  the  disk's  fibers.  The  polariscope  was  equipped  with  the 
standard  devices  for  measuring  isocline  angle,  i.e.  a  pointer  and  gradu¬ 
ated  marks  on  the  circumference  of  the  analyser.  The  disk  could  have 
been  rotated  inplane  by  a  degree  or  two  relative  to  the  zero  degree  dir¬ 
ection  of  the  polariscope.  Related  to  this,  it  was  evident  that  within 
the  disk  the  fibers  were  not  perfectly  straight.  Thus  the  fiber/polari- 
scope  alignment  could  have  been  biased  in  several  ways.  No  effect  was 
assumed  to  be  large.  It  was  simply  assumed  all  effects  added  to  about 


5°. 

In  the  second  numerical  study,  the  exact  optical  isocline  at  each 
mesh  point  was  perturbed  by  adding  to  it  an  angle  that  was  generated  by 
a  Gaussian  random  number  generator.  The  mean  of  these  angles  added  to 
the  exact  angles  was  0°  while  the  standard  deviation  was  5°.  This  ran¬ 
dom  perturbation  represented  possible  experimental  error  in  interpreting 
the  center  of  the  isoclinic  fringe.  The  error  was  assumed  to  result 
from  the  aforementioned  lack  of  good  definition  of  the  isoclinic 
fringes.  All  stresses  for  all  8  fiber/diameter  orientations  were  compu¬ 
ted  and  compared  with  the  finite-difference  computations  using  exact 
values  of  9-j-  and  Nj  of  each  mesh  point. 

As  expected,  a  myriad  of  results  were  generated  from  the  uncertain¬ 
ty  studies.  The  effects  of  the  two  types  of  perturbations  varied  wide¬ 
ly,  depending  on  the  fiber  and  diameter  orientation,  and  the  particular 


stress  component. 


For  example,  the  a 


y 


stress  along  the  diameter  perpen¬ 


dicular  to  the  load,  with  the  fibers  aligned  with  the  load,  was  not  af- 
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fected  by  either  type  of  error.  The  same  was  true  of  the  shear  stress 
for  the  case  with  the  fibers  30°  relative  to  the  load  and  the  diameter 
perpendicular  to  the  fibers.  On  the  other  hand,  there  were  cases  where 
one  type  of  error  had  an  effect  but  the  other  type  did  not.  For  exam¬ 
ple,  the  a  stresses  along  a  diameter  perpendicular  to  the  load  and  with 
the  load  perpendicular  to  the  fibers  were  more  sensitive  to  random  er¬ 
rors  than  to  the  constant  error.  Finally,  some  stresses  were  seriously 
affected  by  either  disturbance  to  the  isocline.  The  shear  stresses 
along  a  diameter  perpendicular  to  the  load  and  with  the  fibers  30  rela¬ 
tive  to  the  load  were  very  sensitive  to  isocline  errors. 

Tables  4  and  5  rank  the  various  situations  as  to  how  they  were  af¬ 
fected  by  the  errors.  The  top  of  the  column  corresponds  to  conditions 
least  affected.  There  is  no  significance  to  the  fact  that  is  the 
left  column,  the  middle  column  and  is  the  right  column.  That  is, 
it  cannot  be  said  that  the  least  affected  stress  was  less  affected 
than  the  least  affected  stress.  It  is  possible  to  say,  however,  that 
for  a  specific  fiber  orientation  and  diameter,  the  larger  of  the  three 
stresses  was  generally  least  sensitive  to  errors.  In  the  tables  the 
figure  number  which  illustrates  the  fiber  orietation  relative  to  the 
load,  and  the  particular  diameter,  is  indicated  in  parenthesis.  Obvi¬ 
ously  these  figures  do  not  represent  the  perturbed  stress  computations. 

The  major  issue,  however,  was;  Could  the  discrepancies  in  the 
stresses  between  theory  and  experiment  observed  in  the  various  situa¬ 
tions  for  the  actual  disk  be  explained  by  the  presence  of  either  a  con¬ 
stant  or  random  error?  Figure  111  shows  the  predicted  a  stresses  along 

A 

the  diameter  perpendicular  to  the  fibers  for  the  case  when  the  fibers 
were  30°  relative  to  the  load.  The  figure  shows  the  values  predicted 
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Table  4 


Ranking  of  Effects  of  Constant  Isocline 
Error  on  Stress  Predictions 
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Table  5 


Rank! ng 

of  Effects 

of  Random  Isocline 

Error  on  Stress 

Predictions 
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using  the  exact  values  of  the  isocline  and  predictions  based  on  the  two 
types  of  disturbances  to  the  isoclineSo  This  figure  should  be  compared 
to  f1g„  9io  The  case  of  the  constant  error  in  fig.  Ill  looks  somewhat 
like  the  mismatches  in  fig,  91,  Figure  112  shov/s  the  predicted 
stresses  along  a  diameter  perpendi cul ar  to  the  fibers  for  the  case  of 
the  fibers  making  a  45°  degree  angle  relative  to  the  load.  This  figure 
should  be  compared  with  fig,  96, 

From  the  figures  it  is  clear  a  constant  error  was  much  more  likely 
than  random  errors  to  explain  the  discrepancies  observed  between  theo¬ 
retical  and  experimentally  measured  stresses.  However,  there  was  a  sig¬ 
nificant  difference  between  the  discrepancies  observed  in  the  actual 
disk  and  the  discrepancies  generated  in  the  numerical  error  studies.  In 
both  fig.  111  and  112  it  is  seen  that  the  constant  isocline  errors 
seemed  to  affect  stresses  near  the  center  of  the  disk  more  than  they  af¬ 
fected  the  stresses  near  the  disk’s  edge.  For  example,  in  fig.  112,  the 
triangles  and  circles  coincide  up  to  about  x/D  =  0.15.  In  fig.  Ill 
these  symbols  coincide  up  to  x/D  =  0.2.  After  those  locations,  the  two 
predictions  diverge  as  the  center  of  the  disk  is  approached.  In  con¬ 
trast,  fig.  96,  for  example,  shows  a  constant  discrepancy  along  the 
v-.'hole  diameter,  starting  at  the  edge  of  the  disk.  This  would  seem  to 
indicate  a  constant  error  of  a  type  different  than  a  constant  isocline 
error.  As  vms  felt  all  along,  and  as  these  isocline  error  studies  con¬ 
firmed,  errors  in  the  level  of  residual  birefringence  were  more  likely 
to  account  for  discrepancies  than  other  types  of  errors. 

For  completeness,  similar  scatter  studies  were  conducted  for  the 
isotropic  disk  stresses.  Recall,  the  experimentally  determined  stresses 
for  the  isotropic  disk  v/ere  quite 


close  to  the  predictions  of  the 


a„/(|P|/Dt) 


Fig.  112  Effect  of  isocline  errors  on  predicted  a  stresses  along 

diameter  perpendicular  to  fibers,  fibers'^45°  relative  to  load. 
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elasticity  solutions^  Figure  109  illustrates  the  stress  which  deviatea 
the  most  from  predictionSo  Figure  113  shov/s  the  stresses  computed  for 
that  same  situation  but  with  the  two  types  of  errors  incorporated  into 
the  isoclinOo  As  can  be  seenj  both  types  of  errors  affect  this  stress 
calculation  to  a  significant  degree,.  Comparing  figs.  109  and  113s  It 
appears  something  other  than  random  errors  caused  deviation  in  fig.  109. 
However  since  neither  type  of  error  caused  the  constant  shift  along  the 
entire  diameter,  as  for  example  in  fig.  89,  the  isotropic  disk  was  not 
studied  further. 

The  results  of  the  error  study  indicated  that  while  isocline  errors 
can  lead  to  serious  errors  in  the  prediction  of  stresses,  such  errors 
were  not  likely  the  source  of  the  errors  observed  in  the  experiments 
with  the  orthotropic  disk.  The  errors  observed  had  very  few  of  the 
characteristics  of  random  errors  and  they  did  not  seem  to  behave  exactly 
like  a  constant  error.  These  errors  were  thus  discounted  as  a  major 
source  of  discrepancy. 
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0.3 


no  error 


/  5°  constant  error 


i  5°  random  error 


O 


disk  diameter  =  D 
disk  thickness  =  t 


Fig.  113  Effect  of  isocline  errors  on  predicted  t  stresses 
along  30°  diameter,  isotropic  disk. 
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SUMMARY  AND  CONCLUDING  COMMENTS 


A  considerable  amount  of  information,  both  in  the  form  of  the  deri¬ 
vation  of  equations  and  in  the  form  of  data,  has  been  presented.  In  the 
second  section  equations  governing  the  photoelastic  phenomenon  in  ortho¬ 
tropic  materials  were  derived.  In  the  third  section,  calibration  of  the 
material  was  discussed  and  calibration  data  were  presented.  The  fourth 
and  fifth  sections  dealt  with  testing  and  checking  the  orthotropic  pho¬ 
toelastic  theory  in  the  context  of  a  unidirectional  disk  in  diametral 
compression.  The  just-presented  section  discussed  sources  of  errors  in 
photoelastic  materials.  As  data  were  presented,  they  were  discussed  and 
put  into  context  with  other  data.  The  purpose  of  this  section  is  to 
summarize  the  important  results  presented. 

The  theory  of  orthotropic  photoelasticity  was  developed  in  the  sec¬ 
ond  section  so  that:  (1)  the  effects  of  residual  birefringence  could  be 
accounted  for  in  the  photoelastic  phenomenon,  and  (2),  the  key  steps 
leading  to  specific  operational  equations,  e.g.  stress  prediction  equa¬ 
tions,  were  shown.  The  experiments  for  calibrating  the  photoelastic  ma¬ 
terial  revealed  that  uniaxial  tests  are  not  the  best  way  to  calibrate 
the  material.  The  four-point  bend  specimen  was  found  to  be  superior  to 
the  uniaxial  specimen  for  calibration  purposes.  The  reasons  were  as 
follows: 

(1)  The  stress  gradients  in  the  bending  specimen  will  lead  to 
sharper  fringes  and  thus  better  fringe  definition. 

(2)  A  single  fixture,  specimen,  and  test  can  be  used  to  load  the 
material  in  both  tension  and  compression. 
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(3)  The  optical  neutral  axis  is  a  convenient  way  to  assess  residual 
bi refringence. 

Unfortunately,  uniaxial  tests  must  be  conducted  for  determining  elastic 

properties  and  so  the  four-point  bend  tests  constitute  an  extra  set  of 
tests. 

The  comparisons  of  the  computer-generated  and  the  observed  fringe 
patterns  revealed  that  on-a-whole,  the  notion  of  orthotropic  photoelas¬ 
ticity,  as  presented,  is  valid.  As  with  the  uniaxial  tests,  sharpness 
of  fringes  was  shown  to  be  a  concern.  Also,  it  was  shown  that  the  level 
of  residual  birefringence  affected  the  correlation  between  the  predicted 
and  observed  fringe  patterns.  By  adjusting  the  assumed  level  of  residu¬ 
al  birefringence  the  predicted  and  observed  fringe  pattern  could  be  made 
to  coincide  for  a  particular  case.  Unfortunately,  the  adjustment  was  a 
function  of  the  load  orientation  relative  to  the  fibers.  The  computer¬ 
generated  fringes  also  showed  that  in  the  presence  of  residual  birefrin¬ 
gence,  the  sign  of  the  applied  load  can  strongly  affect  the  fringe  pat¬ 
tern. 

A  detailed  examination  of  the  predicted  and  observed  photoelastic 
behavior  was  provided  by  studying  the  isochromatic  and  isoclinic  fringes 
along  specific  disk  diameters.  This  type  of  point-by-point  analysis  was 
necessary  for  accurate  determination  of  the  stresses.  The  fringe  pre¬ 
diction  study  was  contrasted  by  conducting  similar  studies  on  isotropic 
disks.  From  the  fringe  prediction  studies,  it  can  be  concluded  that: 
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1)  The  material  seemed  to  exhibit  behavior  associated  with  lower 
rather  than  higher  levels  of  residual  bi refringenceo  Specific¬ 
ally,  a  value  between  0  and  0o07  for  M[^  would  have  correlated 
well  with  80%  of  the  cases o 

2)  Compared  to  isotropic  photoelasticity,  orthotropic  photoelas- 
ticity,  as  hypothesised  here,  does  not  seem  as  accurateo  This 
may  be  related  to  the  lack  of  good  fringe  resolution  or  more 
basically,  to  the  assumptions  of  the  theory  itselfo  Three  as¬ 
sumptions  v^hich  could  be  suspect  are:  (a)  validity  of  through- 
the-thickness  averaging,  (b)  validity  of  the  stress-optic  law 
being  orthotropic,  and  (c),  validity  of  the  superposition  of 
residual  and  mechanically-induced  birefringent  effectSo  To 
challenge  these  assumptions  requires  basic  research  beyond  the 
scope  of  this  study^ 

The  section  on  stress  prediction  reinforced  the  idea  that  residual 
birefringence  in  the  material  was  less  than  N|^  =  0ol5o  Specifically  a 
value  between  0„0  and  0o07  appeared  to  be  correct o  It  is  very  important 
to  note,  however,  in  no  computation  nor  experimental  measurement  v^as 
there  overwhelming  support  for  the  assumption  Nj^  =  0.,  Some  value  of  Mf^ 
could  be  chosen  so  that  correlation  was  good„  Unfortunately,  the  level 
of  Nj^  needed  for  good  correlation  varied  from  one  situation  to  the 
otheTo  Hov/ever,  it  can  be  concluded  that,  in  general,  stresses  can  be 
predicted  with  confidence  using  orthotropic  photoelastic  material o  Un¬ 
fortunately  in  stress  prediction,  and  in  every  other  phase  of  the  corre¬ 
lation  betv/een  theory  and  experiment,  when  the  fibers  were  perpendicular 
to  the  load  the  correlation  between  theory  and  experiment  v;as  pooTo 


Finally,  studies  examining  the  effects  of  isocline  errors  showed 
that  while  isocline  errors  can  lead  to  serious  consequences,  they  were 
not  likely  to  have  caused  the  observed  discrepancies  between  the 
experiments  and  the  theoretical  predictions. 
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APPENDIX  A 


The  form  of  the  dielectric  tensor 


The  dielectric  tensor  can  be  expressed  as 


(A-1) 


where  y,,,  y,,,  and  Xyw  are  components  of  the  susceptibility  tensor. 
Numerically,  the  components  of  the  susceptibility  tensor  are  small 
compared  to  unityo  The  quantity  e^^g  is  defined  as 


®ave 
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(A-2) 


Substituting  eq,  A-1  it  is  seen  that 


'ave 


1  + 


Xx  Xy 


(A-3) 


and  thus  e^yg  is  of  order  unityo 
The  quantity  Ae  is  defined  as 


Ae  =y  ("A— ®-^)2  +  o  (A-4) 


Using  eq,  A-l,  Ae  becomes 
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Ae  = 


(A-5) 


The  quantity  Ae  is  small  compared  to  unity  and  in  particular,  it  is 
small  compared  to  egye*  quantity 


(A-6) 


can  be  approximated,  using  the  binomial  expansion,  as 


4^^  ®ave’  ■ 

Similarly, 
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APPENDIX  B 


Elasticity  Solution  for  Stresses  in  an  Orthotropic  Disk 
in  Diametral  Compression 


The  stresses  in  an  orthotropic  disk  in  diametral  compression  can  be 
represented  by  a  state  of  plane  stress.  If  the  x-axis  is  aligned  with 
the  stiff  (fiber)  direction  and  the  y-axis  is  aligned  with  the  soft  (ma¬ 
trix)  direction,  then  the  z-axis  is  perpendicular  to  the  plane  of  the 

disk.  The  nonzero  stresses  in  the  disk  are  a„,  a,,  and  t  (In  the 

X  y  xy 

discussion  here,  x  is  used  to  denote  fiber  direction.  In  the  text  this 
direction  has  been  denoted  as  IM,  Likewise  the  y  direction  here  is 
denoted  2M  in  the  text.  Using  x  and  y  instead  of  IM  and  2M  is  done 
simply  so  the  subscripts  on  a,  e,  y,  x,  etc,  are  more  familiar.) 

The  three-dimension  equilibrium  equations,  in  rectangular  cartesian 
coordinates,  reduce  to 


and 


'  +  =  0, 

(8-1) 

dx  5y 

da. 

=  0, 

5x  5y 

(8-2) 

Hooke's  Law  becomes 


(8-3) 


The  only  nontrivial  compatibility  equation  is 


d^Y 
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(B-4) 
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Typically  in  plane  stress  problems,  a  stress  function  F(x,y)  is  defined 
such  that 


6^F 


d^F 

‘'y  ‘  ax2 


(B-5a) 


(B-5b) 


a^F 

’’'xy  “  ■  axay  * 


(B-5c) 


With  eq.  B-5,  equilibrium  is  automatically  satisfied.  Using  Hooke's 
Law,  eq.  B-3,  and  eq.  B-5,  the  compatibility  equation  becomes: 
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(B-6) 


Equation  B-6  admits  a  solution  of  the  form 

F  =  F(x  +  py)  ,  (8-7) 

where  p  is  a  constant. 

Substituting  eq.  B-7  into  eq.  B-6  results  in  the  so-called 
characteristic  equation, 

ailli'*  +  355)^^  ^22 


The  solution  given  by  eq.  B-7  is  only  a  solution  if  p  assumes  spe¬ 
cific  values,  namely  one  of  the  4  roots  of  eq.  B-8.  The  roots  of  eq. 
B-8  are  generally  complex,  the  roots  being  denoted  pj^,  P2,  P3  and 
The  complete  solution  of  eq.  B-6  is  thus 


F  =  F^(x  +  p^y)  +  F2(x  +  P2y)  +  F3(x  +  P3y)  +  F4(x  +  p^y)  .  (B-9) 
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The  roots  are  in  pairs,  two  roots  being  complex  conjugates  of  the  other 
two,  i.e. 

=  ^2  and  =  pp  .  (B-10  a,b) 

It  is  convenient,  and  the  basis  of  the  complex  variable  approach, 
to  define  4  complex  variables  and  as  follows: 


Zi  =  X  +  p^y  ,  Zp  =  X  +  ppY  ,  z^  =  X  +  p^y  and  z^  =  x  +  p^y  . 

(B-11  a,b,c,d) 

Obviously  by  eq.  B-10, 

z^  =  Zp  and  Zp  =  z^  .  (B-12  a,b) 

By  eqs„  B-11  and  B-12,  eq,  8-9  can  be  written  as 
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Usually  the  derivatives  on  the  right  hand  sides  of  eq.  B-15  are 


dF  (z  ) 

<t>i(2i)  =  -4-—  (B-16a) 

dFJz  ) 

<^2(Z2)  =  -4z-“  •  (B-16b) 

=  $;(zj  +  <t.'(zJ  (B-17a) 

5x2 


5^F  2  2 

”  ^  ^  1  ^  ^2  ^2  ^  ^  p  ^ 

5y2 

(B-17b) 

s^F 

(B-17c) 

where 

d<t..  (z. ) 

=  --Jz-J- 

(B-18a) 

and 

d(|)p(Zp) 

=  “'dz^^  * 

(3-18b) 

With  this  notation 

2  2 

-  2RG  [  (|)|  (  )  +  j,i^  (j)^  (  Z  r>  )  ]  , 

(B-19a) 

^y 

=  2Re[(j)|(zj^)  +  * 

(B-19b) 

and 

"^xy 

=  -  2Re[p,j^(i)|^(Z]^)  + 

(B-19c) 

The  complex  variable  approach  to  these  plane  stress  problems  is 
convenient  because  (1),  since  it  is  assumed  F  is  an  analytic  function  at 
all  but  a  few  points  in  the  domain,  the  functional  form  of  F  is  known, 
and  (2),  the  boundary  conditions  can  be  related  to  F  through  rather  sim¬ 
ple  integrals.  Concerning  the  latter,  let  and  represent  the  x- 


redefined,  i.e. 


then 
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and  y-direction  surface  tractions,  respectively,  on  the  boundary  of  the 
body,  in  this  case  the  disk. 

Assume  that  on  the  disk  the  normal  surface  traction,  a^,  and  the 

tangential  surface  traction,  are  specified  functions  of  0.  This 

idea  is  shown  in  fig.  B-1.  The  quantities  a-(e)  and  t  (0)  cannot  be 

f  r  0 

completely  arbitrary  in  that  the  disk  must  be  in  translational  and 
rotational  equilibrium  as  a  rigid  body.  With  denoting  the  x- 
component  of  the  surface  traction  and  denoting  the  y-component  of 
surface  traction,  then 

Xn  =  X^(0)  =  (ap(e)  COS0  -  t^q(9)  sin(0))  (B-20a) 

Yn  =  Yg(0)  =  (0^(0)  sine  +  TgQ(9)  coso)  .  {B-20b) 

The  boundary  conditions  on  F  become 

0 

2Re[(^^(z;^)  +  (j)p(zp)  ]  =  -  /  Y  (0)Rd0  +  CT  (B-21a) 

'  '  0 

e 

2Re[^I^(t)^(z^)  +  P2<I'2^^2)]  =  •*■  /  Xg(0)Rd0  +  03,  (8-21b) 

0 

where  and  are  evaluated  at  z  =  Re^®,  i.e.  the  boundary  of  the 
disk. 

Lekhnitskii*  states  the  solutions  to  and 

geometry  are: 

00 

(^l(zi)  =  +  A^z^  +  ^  A^P^  (z;^)  (8-22a) 

m=2  m 

CO 

4.2(22)  =  +  B^zp  +  I  B^P2  (Z2)  .  (B-22b) 

m=2  m 


*Lekhnitskii,  S.  G.,  Anisotropic  Plates,  translated  from  2nd  ed.  by  S. 
W.  Tsai  and  T.  Cheron,  Gordon  a'nd  "Breach,  Science  Publishers,  New  York 
1968,  eq.  33.6. 
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B-1  Surface  tract 


The  Pim(zi)  ^nd  P2m('^2^  ^Iven  by 


(B-23a) 


*"201^^2' 


-1 


'(l-iP2)" 


+  \ 


z^-  R^(l+n^))" 


(8-23b) 


That  these  are  solutions  to  and  '’s  not  obvious  to  even  the 

serious  observer  of  mathematics,  let  alone  the  casual  observeco  These 
solutions  to  and  ^  other  solutions  for  other  domain 

geometries,  is  what  has  made  a  mark  for  Letkhnitskii  and  others  working 
in  this  field. 

Since  (!)j^(Z]^)  and  ipp^Zp)  are  expressed  in  terms  of  a  series,  to 
match  the  boundary  conditions  given  by  eq.  B-21,  Xp  and  of  eq.  B-20 
need  to  be  expressed  as  a  series,  specifically  a  Fourier  series  in  9. 
This  step  is  universal  in  the  complex  variable  approach  to  elasticity, 
i.e.  the  expanding  of  the  known  boundary  data  into  a  Fourier  series. 

Assuming  a^{B)  and  XpgCG)  are  known,  the  boundary  integrals  on  the 
right  hand  side  of  eq,  B-21  can  be  represented  by 

-  Y^(0)Rde  +  (B-24a) 

/  X^(e)Rd0  +  C2  =  Bo  +  I  ’  (B-24b) 

0  m=l 

where 
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a  =  e 


(B-25) 


The  quantities  a^,  Pn,  are  known  coefficients  which  can  be 

determined  by  typical  Fourier  analysis  once  the  functional  forms  of 
X^(e)  and  Y^(0)  are  specified.  Substituting  eq.  B-24  into  the  right 
hand  side  of  eq.  B-21,  and  eqs.  8-22  and  B-23  into  the  left-hand  side  of 
eq.  B-21,  and  equating  coefficients  of  like  powers  of  a  results  in  4 
coupled  linear  algebraic  equations  for  the  A,^  and  B^  of  eq.  B-22.  These 

equations  are: 


m 


-  3. 


m 


+  A„ 


Bn,  = 


(8-26) 


B2^2Bm  ^2^m  "  ^rn 


m  2^3)4....  o 


The  quantities  t^  and  t2  are  given  by 


1  +  i  1  +  i  ij,^ 

1 1  ~  5  and  to  —  1  :r~ 

1  1-1  li^  i-i|j,2 


(B-27) 


For  the  case  m  =  1,  matching  coefficients  in  leads  to 


“1  “1 


Ai  +  +  A^  +  Bj  =  ^ 


(B-28a) 
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(8-28b) 


Equations  828a-c  represent  3  equations  for  4  unknowns,  Ap  8^ 
and  An  indeterminant  solution  is  not  a  problem  if  the  situation  is 
examined  in  detail.  Consider  only  the  linear  terms  in  the  expressions 


for  !})|(zj^)  and  ^2^12)  in  eq,  8-22.  Substituting  only  these 

into  the  expressions  for  stress,  eq.  8-19  results  in 

linear  terms 

Ox  ^  2Re[u^Aj  +  |i^B^  ] 

(B-29a) 

ay  =  2Re[Ai  +  8^  ] 

(B-29b) 

=  -2Re[|i]^Aj  +  « 

(B-29c) 

Eliminating  the  taking  of  the  real  part,  these  three  expressions  can  be 

written  as 

Ox  ■  ^  (B-30a) 

Oy  "  [A|  +  A^  +  +8^] 

(8-30b) 

''xy  ^  ^  ^^2^1  '*■  ^  • 

(B-30c) 

The  right  hand  side  of  eqs.  B-30a,b,c  are  related  directly  to  eqs.  8- 

28c, a, b.  Thus,  the  leading  terms  in  the  expressions  for  stress,  eq.  B- 

30a, b,c,  are 
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Pi  -  Pi 

"^x  “  Ri  ’ 

(B-31a) 

«1  «1 

Oy  -  R 

(B-31b) 

Pi  +  Pi 
\y  -  -  '  “R 

(B-31c) 

Since  only  the  derivatives  of  and  <^>2(22'^  required  for  the 

computation  of  stress,  and  in  eq.  B-22  do  not  enter  into  the  calc¬ 
ulations.  Thus  eqs.  B-31a-c  and  the  solution  to  Eq.  B-26  for  and  B^^, 
m  >  2,  give  the  necessary  constants  to  completely  determine  the  state  of 
stress  in  the  disk.  Evaluating  (})j^(zj^)  and  (t)2(zp)  at  a  specific  x  and  y 
(or  r  and  9)  gives  the  stresses  at  a  specific  location  in  the  disk. 

For  the  diametral  concentrated  load,  the  constants  and 
m  >  0,  are  easily  determined.  As  shown  in  fig.  8-2,  consider  the 
concentrated  load  P  at  an  angle  a  relative  to  the  x-axis.  For  this 
loading 

X^(e)Rde  =  [P6(0  -  a)  +  P6(e  -  a  -  tc)]  coseRde  (B-32a) 
Y^(e)Rde  =  [P6(e  -  a)  +  P6(9  -  a  -  Ti)  ]  sinGRde,  (B-32b) 

0  <  9  <  2-n; 

where  6(  )  is  the 'Dirac  delta  and  +P  represents  diametral  tension.  The 
integrals  on  the  right  hand  side  of  eqs.  B-21  are  1/2-cycles  of  square 
waves,  i .e. 
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/  X  (9)Rde  =  P[U(0  -  a)  cosa  +  U(0  -  a  -  it)  cos(a  -n:)  ]R 
0 

(B-33a) 


0 

/  Y  (0)Rd0  =  P[U(0  -  a)  sina  +  U(0  -  a  -  n)  sin(a  -  it)  ]R  , 

0 

(8-33b) 

U(  )  denoting  the  unit  step  function.  Since 


cos(a  -  Tc)  =  -cosa  and  =  -sina, 

(8-34a,b) 

eq.  B-33  becomes 


0 

/  X^(0)Rd0  =  Pcosa[U(0-a)  -  U(  0  -  a  -  -re)  ]R  (8-35a) 


0 

/  Y^(0)Rd0  =  Psina[U(0  -  a)  -  U(0  -  a  -  Tt)]R  .  (8-35b) 

0 

These  integrals  are  shown  graphically  in  fig.  8-3.  Using  the  usual 
ideas  associated  with  Fourier  analysis,  the  two  integrals  can  easily  be 
expanded  to  a  Fourier  series  of  the  form  eq.  8-24.  The  results  are: 

(B-36a.b) 

(8-36c) 

(B-36d) 

(8-36e,f) 
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For  m  odd; 


Psina  „  _  Pcosa 

«o  =  -  ~2 "  Po  -  • 


“m  "  {sin(ma)  +  i  cos(ma)} 


^m  Tt 


[sin(ma)  +  i  cos(ma) } 


for  m  even; 


Si  ^  Pm  * 


Obviously  only  a  finite  number  of  terms  are  actually  used  in  eqs. 
B-22  and  B-24.  Concentrated  loads  pose  the  most  difficult  convergence 
problems  for  elasticity  solutions,  particularly  if  infinite  series  are 
truncated.  In  reality,  any  finite  number  of  terms  results  in  a  conver¬ 
gence  problem  right  at  the  point  of  application  of  the  load.  So  the 
question  is,  at  points  of  interest  removed  from  the  point  of  application 
of  the  load,  what  are  the  convergence  characteristics?  Numerical  stud¬ 
ies,  based  on  a  FORTRAN  progam  which  coded  eqs.  B-19,  B-22,  B-23,  B-26, 
B-27,  B-31,  and  B-36,  indicated  that  toward  the  center  of  the  disk,  10 
terms  in  the  series  gave  answers  nearly  identical  to  the  series  with  20 
terms.  Twenty  terms  gave  answers  nearly  identical  to  a  40-term  series 
over  85%  of  the  disk.  A  greater  number  of  terms  caused  underflow  in  the 
computer  being  used. 

As  a  check,  a  finite  element  program,  using  216  constant-strain  el¬ 
ements  in  a  complete  disk,  was  used  to  compute  the  stresses.  The  ele¬ 
ment  arrangement  is  shown  in  Fig.  B-4.  There  was  no  particular  reason 
for  using  216  elements  other  than  the  ability  to  be  able  to  compute 
stresses  at  a  large  number  of  points.  The  particular  finite  element 
program  being  used  computed  stresses  at  the  centroids  of  the  elements 
and  thus  many  elements  were  used.  To  compare  the  stresses  generated  by 
the  complex  variable  solution  with  the  stresses  generated  by  the  finite 
element  solution,  the  following  scheme  as  adopted.  For  a  given  load 
direction  relative  to  the  fibers,  a,  the  stresses  were  compared  along 
the  loaded  diameter,  along  a  diameter  perpendicular  to  the  loaded  diame¬ 
ter,  along  a  diameter  in  the  fiber  direction,  and  along  a  diameter  in 
the  matrix  direction.  Such  a  comparison  would  indicate  how  well  the 
stresses  were  being  computed  in  a  highly  loaded  region,  i.e.  moving  from 
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the  center  of  the  disk  toward  the  point  load  along  a  diameter,  how  well 
the  stresses  were  computed  in  a  lightly  loaded  region,  i.e.  moving  from 
the  center  of  the  disk  toward  a  free  edge  along  a  diameter  perpendicular 
to  the  load,  and  whether  the  stress  diffusion  in  the  stiff  and  soft  dir¬ 
ections  was  properly  represented.  The  elastic  properties  of  the  cases 
studied  were  as  follows: 

Ej  =  28.8  GPa  (4.18  x  10^  psi) 

=  9.40  GPa  (1.36  x  10^  psi) 

G^2  =  3.17  GPa  (0.460  x  10^  psi) 

—  0.3  . 

These  values  led  to  values  for  \i  (refer  to  eq.  B-10)  of: 

li]^  =  2.83  i  ,  P2  “  0.618  i  .  (8-37a,b) 

In  the  two  computational  schemes  the  disk  diameter,  D,  was  3.00  in. 
while  its  thickness,  t,  was  unity.  The  comparisons  are  presented  in 
nondimensional  form,  the  stresses  being  divided  by  P/Dt,  where  P  is  the 
applied  diametral  load.  The  stresses  are  presented  as  a  function  of 
their  x  and  y  locations.  The  x  and  y  coordinates  are  nondimensional i zed 
by  dividing  them  by  the  disk  radius.  With  the  particular  finite  ele¬ 
ment,  the  particular  mesh,  and  the  way  a  point  load  was  applied  to  the 
model,  stresses  were  not  computed  exactly  on  the  loaded  diameter,  nor  on 
the  diameter  perpendicular  to  the  loaded  diameter,  etc.  For  both  the 
finite  element  and  complex  variable  solutions,  the  stresses  were  compu¬ 
ted  at  the  centroids  of  the  elements  and  the  element  centroids  were  not 
exactly  on  the  loaded  diameter. 
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Tables  B-1  through  B-4  indicate  the  comparison  for  4  orientations 
of  the  load  relative  to  the  fiber  direction,  ioe«  a  =  0°,  a  =  30°, 
a  =  60°,  and  a  =  90°.  Twenty  terms  were  used  in  the  complex  variable 
solution. 

From  the  tables  it  is  obvious  there  is  good  agreement  between  the 
finite  element  solution  and  the  complex  variable  solution  over  much  of 
the  disk.  Which  is  the  'correct'  solution  is  not  clear.  The  finite  el¬ 
ement  mesh  should  be  finer  near  the  point  load  to  represent  the  singular 
behavior  of  the  stresses  in  that  region.  The  singular  solution  is 
'built-in'  to  the  complex  variable  solution  but  it  is  not  clear  20  terms 
represents  this  portion  of  the  solution.  The  problem  is  not  itnportant 
because  in  the  actual  experiment  the  disk  material  crazed  badly  near  the 
point  where  the  load  was  applied.  Fringes  could  not  be  observed  and 
compared  with  theoretical  predictions.  In  addition,  the  load  obviously 
was  introduced  over  a  finite  region  rather  than  at  a  point,  as  the 
theory  assumes.  Thus  any  fringes  in  the  immediate  region  of  the  load 
would  be  strongly  affected  by  this  lack  of  idealism.  Any  comparison  be¬ 
tween  theory  and  experiment  in  this  region  of  the  disk  would  be  meaning- 
1  ess. 

Based  on  the  comparison  shown  in  Tables  B-1  through  B-4,  the  com¬ 
plex  variable  solution  was  felt  to  be  quite  accurate  and,  compared  to 
finite  element  solutions,  quite  quick  and  inexpensive  to  run  on  the  com¬ 
puter. 
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Table  B-2 

Comparison  of  Finite  Element  (FE)  and  Complex  Variable  (CV) 
Solutions  for  Disk  with  Fibers  30°  to  Load 
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Comparison  of  Finite  Element  (FE)  and  Complex  Variable  (CV) 
Solutions  for  Disk  with  Fibers  60*^  to  Load 
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Determination  of  Stresses  Using  a  Finite-Difference 
Representation  of  Equilibrium  Equations 


In  an  x-y  coordinate  system,  the  plane  stress  equilibrium 


equations. 

in  the  absence  of  body  forces,  are: 

=  0. 

5x  5y 

(C-1) 

,  ."“y  =  0. 

5x  6y 

(C-2) 

The  stress- 

-optic  law  for  orthotropic  photoelastic  materials 

can  be 

written  as 

^  "  C2  C3  "  N^cos(2ftp)  Npcos(20^), 

(C-3) 

(C-4) 

The  constants  Cj  -  Cg  are  calibration  constants  for  the  material,  Nj 
and  0j  are  the  total  isochromatic  fringe  number  and  the  optical  isocline 
angle,  respectively,  and  Nr  and  e,^  are  the  residual  isochromatic  fringe 
number  and  the  residual  optical  isocline  angle.  The  latter  two 
parameters  can  be  thought  of  as  calibration  constants  for  the  material. 

Consider  a  rectangular  array  of  points  in  an  x-y  coordinate  system 
as  in  fig.  C-1.  For  point  i,  j  +  point  a,  the  following 
approximations  can  be  written: 


where 


and 


""x 

=  k  K'^x)i,j+1  -  ^^x^i,jl? 

(C-5) 

^•'xy 

i,j4 

=  1  /!!xy  +  1 

^  i,j+l  i  ,j 

(C-6) 

ZAy  ^^\y  ^i+i  ,j+l  "  ^’’'xy  ^i-1  ,j+l  ^ 

(C-7) 

i.J 

~  Jly  ^^\yVl,j  "  ^\y^i-l.j^* 

(C-8) 

Equation  C-1  now  can  be  approximated,  at  point  a,  as: 


^°x^i,j+l  ~  ^^x^i,j  4Ay  ^^\y^i+l,j  "  ^-^xy  Vl,j 
(^\y^i+l,j+l  ■  (^"^xy  ^i-l,j+l  ^ 


.  + 


(C-9) 


Equation  C-9  is  an  approximation  to  the  equilibrium  condition  in  the  x 

direction.  It  is  applied  at  specific  points  along  the  midrow  of  the 

three-line  mesh,  row  i.  This  expression  averages  shear  stresses. 

1 

For  point  (i  +  following  approximations  may  be 

written : 


where 


and 


i4,J 

Ay  ^^°y^i+l,j  ’ 

(C-10) 

""xy  J 

5x  i 9  j 

_  1  ^\y  , 

(c-11) 

ax  i+l,j 

"  2ax  ^t^^xy  ^i+l,j+l  ■  ^\y^i+lJ-l^ 

(C-12) 

ax  i,j 

"2^ 

(C-13) 

Equation  C-2  can  be  approximated  at  point  b  as 
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(C-14) 


J  ^  4ax  ^i+l  J+1  "  ^i+l  J-1 

’’’  (^'^xy^ijj+l  "  ^^xy^i,j-l^  ” 

Equation  C-14  can  be  applied  at  discrete  points  b  and  b'  along  the 
columns  of  the  mesh,  columns  j-1 , . . . ,j+2 . . «  The  equation  also 

averages  shear  stresses. 

Equations  C-3  and  C-4  can  be  written  as  follows  to  put  them  into 
context : 


.  iilu  +  Ky_!i,l  .  N.  costae,  )  -  N„cos(29„1 
^2  4  'i,j  'i,j 


(C-15) 


^^x^j,j  '  ^'"'xy ^i , j  _  ^i , j 


sin(2e^  )  -  y-  sin(2e  ). 


(C-16) 


(Here  N;^  and  0)^  are  assumed  to  be  the  same  at  each  mesh  point). 

Equations  C-9,  C-14,  C-15,  and  C-16  are  the  four  equations  which 
can  be  used  to  compute  the  stresses  Oy,  and  at  each  point  i , j . 

It  is  obvious  that  x  direction  equilibrium  is  not  enforced  at  the  same 
points  as  y  direction  equilibrium  and  here  y  direction  equilibrium  is 
enforced  at  twice  as  many  points.  There  is  nothing  special  about  this, 
it  is  simply  a  direct  result  of  the  mesh  geometry,  i.e.  long  and  thin. 
The  photoelastic  equations  should  be  applied  at  every  point.  The  equil¬ 
ibrium  equations  can  be  applied  so  there  are  exactly  as  many  equations 
as  there  are  unknown  stresses.  Alternatively,  the  equilibrium  equations 
can  be  applied  frequently  enough  so  that  there  are  more  equations  than 
there  are  unknown  stresses.  Solving  an  overdetermined  set  of  equations 
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will  give  stresses  that  satisfy  all  four  equations  in  a  least-squares 
sense. 

Finally,  of  course,  there  are  boundary  conditions.  At  some  points 
one  or  more  stresses  may  be  known,  e.q.  (<7^)4  5  “  These  are  simply 
more  equations  in  the  equation  count  and  can  be  used  to  give  a  unique 
solution  or  an  overdetermined  solution. 
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